ru/\sA 

Technical Memorandum 86183 


21 Layer Troposphere- Stratosphere 
Climate Model 

^ (NASi-TM-86183) 21 1AYEB 1185-21671: 

TBOBOSPHBEE-STfiiTOSPHIBE CXIKATE MODEL 
- (MSA) 78 p HC AC5/ME A0 1 CSCL 04B 

Onclas 

I JiyS 

D. Rind, R. Suozzo, A. Laois, G. Russell, 
J. Hansen 


December 1984 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 
Institute for Space Studies 
New York, New York 10025 



ABSTRACT 


Our global climate model, (Model II). is extended through the stratosphere 
by increasing the vertical resolution and raising the rigid model top to the 
O.Olmb (75 km) level. The inclusion of a realistic stratosphere is necessary 
for investigation of the climate effects of stratospheric perturbations, such 
as changes of ozone, aerosols or solar ultraviolet Irradlance, as well as for 
studying the effect on the stratosphere of tropospheric climate changes. 

The model is Integrated for 15 months with 8* x 10® resolution and 21 
layers. Results are shown for the full seasonal cycle. The model 
realistically simulates the observed temperature and wind patterns throughout 
the troposphere and stratosphere. This Includes latitudinal gradients, such 
as the reversed temperature gradient in the winter lower stratosphere and the 
summer stratosphere; seasonal variations, such as the change from winter 
westerlies to summer easterlies, and semiannual oscillations in temperature 
and wind in the equatorial upper stratosphere (although the west wind phase is 
weak); and realistic hemispheric differences, such as the colder lower winter 
stratosphere in the Southern Hemisphere, the greater zonal kinetic energy in 
the Southern Hemisphere %rf.nter stratosphere, and the greater eddy kinetic 
energy in the Northern Hemisphere winter stratosphere. The Southern 
Hemisphere has a larger transient to standing eddy kinetic energy ratio in 
winter than the Northern Hemisphere. The simulation develops what appear to 
be Inertial oscillations in the equatorial upper stratosphere of the %7lnter 
hemisphere. 

The model also realistically reproduces the observed wave 1 and wave 2 
forcing at the base of the lower stratosphere, and the long wave amplitudes 
and phases are similar to observations throughout most of the stratosphere. 

In the upper stratosphere the long wave amplitudes are too strong, possibly 
due to the lack of photochemical damping. Energy budget calculations show 
that the long waves are producing an energy cycle while propagating through 
the lower and mid stratosphere of the winter hemisphere. Eastward moving 
tropical waves display many of the properties expected of Kelvin waves 
throughout the stratosphere and appear to be associated with the semiannual 
wind oscillation. Horizontal transports of sensible heat and westerly 
momentum are realistic throughout the troposphere and stratosphere, except for 
excessive values in the upper stratosphere associated with the large wave 
amplitudes. The vertical transport of geopotential energy and the 
Ellassen-Palm (EP) flux show the expected greater penetration in the winter 
hemisphere, and the EP flux is directed into the core of the stratospheric jet 
due to the relatively large gradient of quaslgeostrophlc potential vortlclty 
encountered there. The EP flux divergence is similar to that observed 
throughout the atmosphere. Finally, a time sequence of downward propagating 
minor warmings is shown which displays many of the properties associated with 
actual stratospheric wainnlngs. 

In addition to the excess planetary wave amplitude in the upper 
stratosphere, other model deficiencies Include the Northern Hemisphere lower 
stratospheric temperatures being 5-10®C too cold in winter at high latitudes, 
probably due to the model's failure to allow excess water vapor to condense 
above lOOmb. Also, temperatures at 50-60km altitude near the equator are too 
cold, probably as a result of excessive CO 2 cooling. This latter effect is 



probably also responsible for the temperature gradient being too large at 
those altitudes near the equator, and thus for east winds of large magnitude 
in the tropical mesosphere* Methods of correcting these deficiencies are 
discussed* 



!• Introduction 


It is now recognized that the stratosphere may have a strong Influence on 
lower atmospheric dynamics and climate on all time scales. For example, there 
Is the possible Influence of the top boundary condition on weather forecasts 
(Llndzen et al. , 1968; Spahr et al. , 1982) (thus the possible Influence of the 
stratosphere on actual weather); the effect of the stratospheric refractive 
Index on wave propagation from the troposphere, providing a mechanism for 
solar cycle effects on weather and climate (Hines, 1974); on a longer time 
scale, the changes that the ozone amount or vertical distribution would have 
on climate (Ramanathan and Dickinson, 1979). These effects are often 
Interactive: for example. If changes in ozone were to change the tropospheric 

temperature, the water vapor coming Into the stratosphere might also change, 
further altering the ozone. 

Similarly, changes Initiated In the troposphere may have strong Impacts 
on stratospheric structure and composition. Short time-scale events such as 
blocking have been Implicated In the generation of stratospheric warmings 
(Craig and Bering, 1959). Warming of the troposphere due to Increased carbon 
dioxide could plausibly alter the water vapor content of the stratosphere, 
with the subsequent chemical alterations affecting numerous species (this Is 
in addition to the direct effect of added CO 2 cooling In the stratosphere). 
These types of events are In addition to the well-known Influence of the 
troposphere on the dynamics of the winter stratosphere. 

In order to evaluate and understand these effects in a numerical model it 
is necessary that the model realistically simulate both tropospheric climate 
and the stratosphere. This paper describes one such model; it is an extension 
of the climate model described by Hansen et al. (1983 - hereafter referred to 
as I). There have been only two other primitive equation troposphere- 
stratosphere models published In the literature (Fels et al. , 1980; Hunt, 

1981). This model differs from both In that it Includes full seasonal and 
diurnal cycles, and does not Include explicit diffusion. The model contains 
all the detailed calculations and parameterlzatlons necessary for climate 
simulation, as given In I. The differences between this model and the 9 layer 
"model II" discussed In I mainly concern the numerical techniques or 
parameterlzatlons In the stratosphere; these are discussed In Section 2. 

Selected output from the model is shown and compared with observations In 
Section 3. There are several purposes behind this comparison. The most 
obvious is to gage the acceptability of the model In simulating the current 
troposphere and stratosphere. To this end the comparison with observations Is 
extensive, similar to that which was done In I for the tropospheric simulation. 
This has not previously been published for a full troposphere-stratosphere 
model, and is made possible by the continued accumulation of observations 
throughout the stratosphere and climatologies in the troposphere. In the 
course of the comparison opportunities arise to compare different sets of 
observations, and the degree of agreement among them is noted. Problems and 
choices which arise In the course of model development are discussed, again 
with reference to observations. Specific aspects of stratospheric dynamics are 
explored through the model simulations, including semiannual oscillations, 
apparent Kelvin wave variations, the energy cycle at various levels, the 
relative importance of transient versus stationary kinetic energy, and the 
relevance of stationary conservative wave theory and quasl-geostrophlc 



diagnostics In a time-varying primitive equation model. Finally, a series of 
minor warmings Is discussed, which developed In the model during January, the 
first time such events have been shown In this type of model. 

Model deficiencies are discussed as they appear from the comparison with 
observations. They are summarized In Section 4, along with a brief discussion 
of needed Improvements. 

2. Model description 

In this section we will describe the model, emphasizing the differences 
between this version and the version published In I. The model structure Is 
discussed In Section 2a, the fundamental equations In 2b, the numerics In 2c, 
radiation In 2d, convection and condensation In 2e, and other model physics In 
2f. 

a. Model structure 


The model Is global In horizontal extent. The horizontal and vertical 
resolutions are variable. The resolution employed here Is 8° x 10° (latitude 
X longitude) with 21 levels In the vertical, as shown In Fig. 1. (As 
altitudes and pressures will be used Interchangeably In the text, this figure 
provides the equivalence). There are approximately three levels In the 
boundary layer (although this varies with the diurnal cycle), an additional 6 
layers up to the tropopause, 8 layers in the stratosphere, and 4 In the 
mesosphere. For radiation calculations there are an additional three layers 
at the top which do not Interact with the dynamics but serve to provide a 
radiative equilibrium boundary at the top of the atmosphere. In the 
troposphere above the boundary layer, the layers are of nearly constant 
pressure thickness, while above the tropopause they are of nearly constant 
height thickness (log pressure). The vertical coordinate In the troposphere 
Is the a coordinate system (Phillips, 1957) so the ground Is a coordinate 
surface, while at altitudes above lOOmb the model uses constant pressure 
coordinates. This Is necessary In the stratosphere to avoid spurious 
oscillations Induced by a coordinates In regions of variable topography. For 
example, meridional winds In the equatorial region near the Andes show large 
several grid point oscillations when a coordinates are used throughout. 

The model structure at each grid point Is shown In detail In I. Each 
grldbox has aiJ{)roprlate fractions of land and ocean, with the ocean 
temperature and sea Ice specified climatologlcally based on monthly mean 
values. Antarctica, Greenland and some Arctic Islands are covered by Ice 
sheets; snow Is calculated by the model. The land has realistic topography. 
At each grid box radiative fluxes, convection and large-scale clouds are 
calculated, as well as sensible and latent heat fluxes between the ocean or 
land and the atmosphere. Temperature and soil moisture are calculated at 
several levels %rlthln the ground. In these respects the model corresponds to 
"Model II" as published in I; see I for more details. 

b. Fundamental equations 

The equations used are those given In I. They represent conservation of 
mass, momentum, energy, and water vapor, along with the equation of state. 
They differ In the stratosphere from those shown in Table 2 of I only due to 
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the use of constant pressure coordinates. The boundary conditions are that 
a = 0 at o = 1 (the surface), and u = 0 at p “ .Olmb (the model dynamic top). 

Two additional points should be ^de about these equations. First, there 
is no explicit diffusion employed, either in the vertical or the horizontal. 
Thus this model differs substantially from other such models (l.e.. Hunt, 

1981; Pels et al. , 1980). Second, there is no equation for ozone 
photochemical feedback, the ozone values being prescribed. These points are 
discussed further in subsequent sections. 

c. Numerics 


As discussed in I, the model uses the B grid of Arakawa (1972), except 
that potential temperature is a prognostic variable, and the Coriolis force 
and metric term are included at the pole grid points. In addition, in this 
paper there is a major change in the advection term in the energy equation. 

The second order scheme used for the advection of heat in 1 has a tendency to 
produce rather noisy patterns (Russell and Lemer, 1981). This deficiency is 
not very noticeable in the troposphere, partly because the temperature 
gradients tend to be small, and partly because convection mixes heat rapidly 
in the vertical, destroying horizontal noise patterns. In the stratosphere 
the problem becomes severe. With the decreased density, wave amplitudes (wave 
1 primarily) become quite large, as do horizontal temperature gradients 
associated with the waves. The stability of the stratosphere largely prevents 
any convection. Thus the noise in the temperature patterns produced by the 
second order scheme becomes Important, and helps produce local wind gradients 
of large magnitude. These in turn cause difficulty in integration, especially 
in a model such as this with no explicit diffusion. 

To alleviate this problem, the second order differencing scheme was 
replaced in the stratosphere (above lOOmb) by the "slopes scheme” (Russell and 
Lerner, 1981), which produces much smoother patterns in a non-dlf fusive 
manner. In one dimension, this is a higher order upstream scheme, which 
represents the distribution of the advected quantity (in this case, potential 
temperature) within the grid box as a mean concentration and a slope. When 
advection takes place in two or three dimensions, the slopes may be advected 
laterally with an upstream scheme, as in Russell and Lerner (1981), or with 
second order differencing. We tried the slopes scheme in all three directions 
throughout the model with the spatial leapfrog splitting technique (SLF) 
described by Russell and Lerner. This resulted in much smoother patterns of 
temperature and thus height and wind, and allowed the model to be run without 
explicit diffusion. The derivation of this schen», its exact formulation, and 
the results obtained using it compared to second and fourth order schemes, are 
given in Russell and Lemer (1981). 

When applied to the troposphere, this scheme altered the available 
potential energy, and thus the eddy energy and transports. It was desired, 
however, for the model to produce a climatology similar to that published in 
I; thus it was necessary to keep the second order scheme in the troposphere. 
The model thus has a different numerical heat advection scheme in the 
horizontal above and below lOOmb, which is also the level at which it switches 
from a to constant pressure coordinates. 



The possibility arose that using a different vertical advection scheme in 
the troposphere and stratosphere would introduce numerical noise at the lOOmb 
level* Tests were run using only the vertical portion of the slopes scheme in 
the troposphere; they showed that this did not significantly change the 
tropospheric energetics or climate from that shown in 1. Thus the vertical 
portion of the slopes scheme is used in the model throughout the atmosphere, 
and provides continuity through the lOOmb level. The model results give no 
evidence of numerical perturbations generated at or near this level. 

The bottom boundary condition in the model allows for zero vertical 
velocity at the surface, appropriate in a coordinates. The top of the model 
presents a well-known difficulty: the use of an artificial lid traps energy 

which ought to propagate through. The result is a region of high velocity 
winds near the top; from the thermal wind relationship, the large shears 
require strong latitudinal temperature gradients. This is accomplished by 
Inducing rising air at high latitudes, establishing an artificial circulation. 
Reflection of waves off the top should also affect wave transports (Kirkwood 
and Derome, 1977), limiting poleward heat transport. 

Models with tops at around 30km and realistic topography have difficulty 
closing off the tropospheric jet (Tenenbaum, 1982). As the top of the model 
Is raised, the region of excessive winds rises with it; when the top is near 
50km, a level which ordinarily has strong «d.nds, it becomes difficult to 
numerically integrate for any length of time. When the model top is raised 
out of the stratosphere, it then becomes difficult to close off the 
stratospheric jet (e.g., Mahlman and Sinclair, 1979). 

This represents, at least to some extent, a numerical modeling 
difficulty. In addition, at mesospheric heights the induced dlabatlc 
circulation from summer to winter pole, theoretically on the order of 1 m s“l, 
when acted upon by the Coriolis torques, should produce mesopherlc zonal 
motions more than an order of magnitude larger than observed. One solution to 
both of these problems is to Introduce friction into the mesosphere below the 
model top (Leovy, 1964; Schoeberl and Strobel, 1978; Holton and Wehrbeln, 
1980). This technique works to alleviate the problem of energy accumulation 
near the artificial lid, as was shown in I, when it was applied to the 30mb 
level. It also can be used to explain why the observed mesospheric winds 
decrease with altitude even in the presumed presence of the dlabatlc 
circulation. The friction is thought to result from breaking gravity waves 
(including tides), the adiabatic temperature changes associated with the waves 
producing an unstable lapse rate when combined with the already steep 
temperature decline with height prevalent in the mesosphere. If these waves 
have a low level source their (pseudo) momentum is significantly less than the 
ambient momentum, and as the waves break they would act as a drag on the flow. 

In reality, this process, and its application, may be ouch more 
complicated. The point of absorption of these waves may be difficult to 
determine since it may Involve critical level absorption. For this and other 
reasons, the parameterization of the process should differ from the simple 
Rayleigh friction often used (Llndzen, 1981). There is even some question as 
to whether other processes, such as Inertial instability (Dunkerton, 1981) 
might not actually be more responsible for drag in equatorial regions. 



In this model It Is assumed that the dominant drag mechanism In the 
mesosphere is that associated with breaking gravity waves. Frictional drag Is 
Incorporated In layers 19-21 (approximately 65-75km); although analysis shows 
that the effects may occur as low as 50km in winter (Lindzen, 1981), 
observations show that the apparent turbulence is largest above 70km (e.g. , 
Rottger, 1980). The form of the drag is similar to that employed at the 
surface (I, equations 44-49), except the variation of the drag coefficient in 
conditions of non-neutral stability has been altered. The stability criterion 
is determined by assuming that a gravity wave with a surface temperature 
perturbation of 0. 1°C has propagated up to these heights with no attenuation 
(the kinetic energy density remains constant). The wave is assumed to have a 
1km vertical wavelength and induces a peak to peak temperature perturbation 
over a 500m vertical distance. This lapse rate is added to the background 
lapse rate. Interpolated to this scale; if the resulting total lapse rate 
exceeds the dry adiabatic, the drag coefficient is maximized. Insofar as the 
density decreases with altitude, the wave amplitude will be highest at the 
highest levels. The instability will also be largest in regions where the 
ambient lapse rate is largest, for constant wave amplitude. This 
parame-terlzation thus is an attempt to mimic the conditions under which 
breaking gravity waves are most likely. 

The decay time associated with this drag in the winter hemisphere 
averages about 2 ± 1 days at 65km, about 1 ± 1/2 day at 70km, and about 
1/2 ± 1/4 day at 75km (varying with wind speed, and thus latitude). These are 
somewhat faster than those employed by Holton and Wehrbeln (1980) in their 
symmetric model (in which the decay times ranged from 5 days to 2 days at 
these altitudes), although the use of a more realistic infrared radiation 
algorithm would require that model to use stronger mechanical dissipation, 
more like the values used here (Wehrbeln and Leovy 1982; Holton, 1983). 

Various values and forms of the drag were tried, including some with longer 
decay times. The winds and temperatures in the stratosphere were only 
slightly altered by the different magnitudes. 

An additional question which arises is what to do with the energy removed 
by the drag. The same question can of course be asked about the energy lost 
to friction at the earth's surface, but that is insignificant compared to the 
energy in the troposphere. In the mesosphere, however, the results are quite 
different. Several runs were made in which the energy removed by the drag was 
put back in the form of heat, to be lost eventually by radiative fluxes. The 
resulting temperature differences at 70"N, where the drag is strongest in 
winter, between these runs and runs in which the energy was simply discarded, 
ranged between 10®-30®C, with the highest values at the greatest altitudes 
(75km). At 70*’N including the energy back into the model produced results 
which were strongly at variance with the (somewhat uncertain) observations in 
the mesosphere. At other latitudes, where the effects were less, it was not 
obvious which formulation compared better with observations. It also appeared 
as if the change in stability resulting from the different formulations 
affected the long wave amplitudes throughout the stratosphere, although this 
could not be proved as the model has not been run sufficiently to determine 
natural standard deviations. The final version of the model discards the 
energy removed by the drag; despite the theoretical arguments in favor of the 
effects of breaking waves, and the observations of (sporadic) turbulence, the 
inclusion of the drag in the form used here must be considered an artifact of 
the numerical modeling procedure to reduce the influence of the model top. We 


thus in effect hypothesize that the waves propagate through the model's 
artificial lid, and are unavailable to influence events below. 

As mentioned in the previous section, the model Includes no explicit 
diffusion, either in the horizontal or in the vertical. Horizontal turbulence 
likely exists at all levels (e.g. , Zimmerman and Keneshea, 1981), and vertical 
turbulent transfers are likely in the mesosphere (Allen et al. , 1981). The 
exclusion of diffusion is based on the desire to present a model in which 
diffusion is not necessary for numerical stabilization, and does not influence 
model simulations in an unwarranted fashion. Subsequent versions of the model 
may incorporate realistic diffusion values if they can be accurately 
determined. 

The time differencing is the same as that in I, except where it is 
modified in conjunction with the slopes scheme (see Russell and Lerner, 1981). 
The smoothing and conservation of Integral properties in this model are the 
same as those in I. 

d. Radiation 


The radiation routines, both short wave and long wave, are essentially 
those used in I, except for the following differences. 

1) For the short wave radiation, the optical thicknesses for C02» H 2 O 
and ©2 were interpolated directly from numerical k-distrlbution tables as is 
.the case for thermal radiation in 1. (Solar radiation in I uses analytic 
formulas that were fitted to the k-dlstributlons). This improved the accuracy 
in those cases where the analytic fit was imperfect, particularly in the upper 
stratosphere. 

2) In the long wave radiation, the tabulated absorption coefficients were 
refitted to obtain a more accurate integration over zenith angle and closer 
agreement with ID model results for stratospheric cooling rates. Both the ID 
and 3D radiation codes use explicit integration over zenith angle; since 
zenith angle effects become more important in the stratosphere, there is some 
need to compensate for the coarser 3D resolution. This change is more 
Important for improving the accuracy in the more rarified regions of the 
atmosphere. 

3) In the long wave radiation, the 3D code merges the spectral and 
probability intervals of the k-dlstrlbutlons compared to using an explicit 
spectral integration scheme with typically 10 probability intervals per 
spectral band in the ID model. We refined the parameterlzations for merging 
the overlapping CO 2 and water vapor absorption to provide closer agreement 
with ID model results for the downward directed thermal fluxes. This change 
basically Improved the details of the surface energy balance, but is of lesser 
importance in the stratosphere. 

One other change involved the way the temperatures are transmitted to the 
radiation scheme. The use of the slopes scheme in the vertical provides a 
first order approximation of the temperatures at the edges of the layers, in 
addition to the normally calculated value for the layer as a ;%diole. As the 
radiation assumes a linear Planck function within a layer, a slope for oT^ was 
least-squares fit to the linear T slope produced by the slopes scheme in the 



dynamics subroutine. This then defined the temperature at the layer edges 
with more resolution than would have been the case otherwise. 

The column average climatological ozone distribution as a function of 
latitude and month was taken from London et al. (1976). The vertical 
distribution is also a function of latitude and month, at altitudes below the 
lOmb level based on the distributions given by Dutsch et al. (1971, 1974) and 
at altitudes above lOmb from London et al. (1977). The resultant distribution 
is generally consistent with the more recent observations (McPeters, 1980; 

Wang et al., 1982), although they do not maintain the hemispheric asymmetry in 
upper stratospheric ozone (e.g. , Maeda and Heath, 1973). The ozone 
distribution has not been altered to offset other radiation errors by tuning 
the equatorial stratospheric temperatures as has been done in other models; 
the philosophy is to present the model results unaltered so as to highlight 
uncertainties in modeling or theory. In addition, as noted previously, the 
ozone does not respond photochemlcally to changes In temperature. Although 
the ozone climatology Includes longitudinal variations, in the model we used 
zonally averaged ozone amounts. 

The water vapor distribution at altitudes of 18km and above was 
initialized at a value of 3 X 10~6 {| 20 /kg air), a value consistent with 

observations (WHO, 1981) although there is no Increase with altitude 
prescribed. Due to the uniform distribution, after 15 months values 
throughout much of the upper stratosphere and mesosphere were close to the 
initial value. Thermal cooling would be slightly different if other values 
were used in the Initialization. 

In the lower stratosphere large-scale advection alters the Initial water 
vapor distribution, due at least partially to unrealistic restrictions 
discussed below. 

e. Convection and condensation 


The only difference between the moist convection scheme outlined in I and 
the version employed here relates to the transition level between a and 
pressure coordinates. Because of the complexity Involved in the scheme, it is 
difficult to calculate the mass flux through this transition; thus in this 
model moist convection is not allowed to pass through the lOOmb level. In 
addition, lar^-scale condensation is also not allowed above lOOmb, again for 
the sake of simplicity. In earlier versions of the model, which were 
completely in a coordinates, moist convection and condensation were allowed to 
occur at all altitudes, and the mass flux through lOOmb, and condensation from 
100-50mb, was small but not zero. The restrictions Imposed here thus produce 
errors in the moisture profile. There is nothing about these schemes that 
inherently prohibit their use in a and pressure coordinates, and these 
restrictions should be removed in future model integrations. 

Dry convection is allowed at all levels, 
f . Other model physics 

The rest of the model physics is the same as that discussed in I. This 
Includes cloud calculations, ground physics, and surface air quantities. The 
model computes clouds at all levels, temperature and moisture at two levels 



within the ground, enow cover, and surface air temperature, wind and moisture 
for flux and diagnostic purposes. The version used In this model specifies 
the sea surface temperature and ocean and land Ice cover (see I for more 
details). 

3. Model Simulations 

The model described above was Integrated for 15 months, starting from 
climatological Initial conditions for December 1. The boundary conditions 
affected by the seasonal cycle (e.g. sun, sea surface temperatures) were 
varied dally, and a diurnal cycle was Included. The results presented here 
will be for the last 12 months of the simulation. The results for February of 
year 2 were similar although not Identical to those for February of year 1, 
implying that the three month adjustment was sufficient for most quantities. 

In the troposphere, even with specified sea surface temperatures, there are 
longer time constants In the system (such as ground moisture and ground 
temperature in the second layer) for which three months are not sufficient. 

As the tropospheric simulations were discussed In detail In I, we will 
concentrate on the results for the stratosphere, except where results In the 
troposphere differ from those In the 9 layer model, or are important for 
understanding the stratospheric simulation and were not previously presented. 

Comparison will be made with observations for a number of different 
quantities. This presents several problems. The observations often differ 
substantially from one another, either due to the use of different techniques 
' or to the use of different years of data. Especially in the Northern 
Hemisphere winter, the Inclusion of years with stratospheric warmings produces 
zonal averages which vary widely from those compiled for years %d.thout major 
warmings. The uncertainty In the meaning of a "typical" winter also may apply 
to the model: the model winter used for this comparison did not feature a 

major warming and there is no way to know what the Interannual variability in 
the model will be without further Integrations. We will thus try to compare 
the model results with those for typical non^armlng winters, although the 
collected climatologies combine both types of winters. 

For many quantities, especially in the Southern Hemisphere, we are 
limited to one or perhaps two years of analysis at most. These comparisons 
can thus not be thought of as definitive* 

a. Temperature 

The model simulations will be shown In a variety of formats, for both 
January and July* 

1) ZONAL MEAN TEMPERATURE VERSUS ALTITUDE 

Fig. 2 displays the results for January and July of the zonally averaged 
temperature as a function of latitude and pressure. Realistic features 
Include temperature decreasing with Increased altitude In the lower 
stratosphere, increasing above to the stratopause, then decreasing in the 
mesosphere. The summer stratosphere and troposphere are warmer than the 
winter ones, while the reverse Is true in the mesosphere. A more detailed 
analysis Is given below. 
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Fig. 2. Latitude-height profile of zonal mean temperature for (a) January 
and (b) July* Dashed lines indicate in-between contour values* 




JANUARY: The tropical tropopause temperature of -74®C compares 

favorably with the value of -76"C reported by Geller et al. (1983) for the 
years 1978-1981. This Is an Important parameter If the region Is acting as a 
cold trap for water vapor entering the stratosphere. It Is also In agreement 
with the reports of Newell et al. (1969), but Is somewhat wanner than the 
-80“C given by Oort and Rasmusson (1971). As pointed out by Newell (1972) the 
variation of this temperature Is probably related to changes In the vertical 
motion associated with the Hadley circulation (whose estimates are also quite 
variable, as noted In I, and which may be subject to large year-to-year 
variations). Variations of several degrees also occur In this region due to 
the quasl-blennlal oscillation and volcanic aerosols (e.g., Angell and 
Korshover, 1964; Quiroz, 1983), as well as the solar cycle (Quiroz, 1979). 

(The model value may also be affected by the vertical layering, as temperature 
changes rapidly with altitude). 

The temperature gradient In the model's lower stratosphere properly 
reverses from the equator to 55“N, but the values north of about 60 “N are too 
cold by 5-10“C. One possible contributing factor Involves the water vapor, 
which as noted In Section 2e was not allowed to condense above lOOmb. Thus In 
January the specific humidity poleward of 50°N was six times too large at 
73mb, three times too large at 34mb, and two times too large at 16mb compared 
with observations of 3-4 ppmv (see the compilation of observations In the WHO 
report, 1981). Using Standard Atmosphere temperatures for a one time-step 
Integration, this water vapor produced cooling rates which ranged from two 
times too large at 73mb, to 15Z excessive at 16mb. A two month experiment was 
conducted to determine the effect the proper water vapor would have on the 
model's temperatures. Starting with the initial conditions for December 1, 
December and the final January were rerun with the water vapor fixed at the 
observed value. The monthly average temperatures In January were 2-3“ warmer 
in the high latitude lower stratosphere, and seven degrees warmer In the 
middle stratosphere, while the cooling rate for the month was still reduced by 
some 33Z at 73mb. It thus seems likely that the excessive water vapor, which 
would have been reduced If condensation had been allowed (as average relative 
humidities were well above 100 %), was greatly If not entirely responsible for 
the cold temperatures In the model simulation. Another potential contributing 
factor, sensible heat transports, will be discussed in Section 3e. 

The excessive water vapor peaks at the 150 and 73mb layers, and then 
decreases to observed values by 7mb. At 30mb high latitude winter 
temperatures are still slightly too cold compared to climatological averages 
(e.g., Geller et al. , 1983) but are actually In good agreement with "cold" 
January observations, such as In January of 1981 (Labltzke and Goretzkl, 1982). 
At the lOmb level the temperature throughout the Northern Hemisphere Is In good 
agreement with a typical "cold" January (Labltzke, 1981), and It remains so 
through the 1.6mb (45km) level (Geller et al. , 1983; Groves', 1971). 

In the region extending from 50-60km at low latitudes, temperatures once 
again are colder than observed, by 15-20“C. As this region Is not strongly 
affected by dynamic transports, the problem was thought to Involve the radiation 
calculations. A comparison was made between the CO 9 cooling rates calculated In 
this model, and those reported by Apruzese et al. (1982). With the radiation 
scheme employed here the cooling rate from CO 2 increased from near 6 “C day~^ at 
45km to 9-10“ day“l at 50-56km for the D.S. Standard Atmos-phere (1962) 
temperature profile. In the several models reported by Apruzese et al. (1982) 



the cooling rates were significantly lower from 50-56km, ranging from 7-8“C 
day”l. While the lack of agreement in Itself would not prove there is excessive 
cooling in this radiation scheme, the fact that the temperatures are so low 
compared with observations makes it likely that the cooling above 45km is 
overestimated. The current k-dlstrlbutlon method was developed for the more 
massive lower atmosphere; work is in progress to improve the scheme for the 
region above 45km. 

The relatively warm high latitude temperatures observed, near 50km at high 
northern latitudes during this month result from warming pulsations which 
occur associated with the large amplitude wave number 1 In the model. This 
will be discussed In Sections 3d and 3f. 

The Southern Hemisphere model values during this month appear to be In 
good agreement with available observations (Newell et al., 1969; Koshelkov, 

1977) although they might be 5“C too cold in the polar upper troposphere. 
However, the successful simulation of the extratroplcal summer stratospheric 
temperatures at 50km, In contrast to the relatively cold equatorial 
temperatures at that altitude, raises the possibility of compensation due to 
excessive ozone heating at Southern polar latitudes. To Investigate this the 
model ozone heating rates were compared with those calculated by several 
modelers for the 50km region. On January 1 the heating rate In the model 
varied from about 9“C day“l at the equator to around 21 ®C day“^ at 90“S. This 
Is In essential agreement with the results of Park and London (1974) from whom 
the ozone distribution above lOmb was derived. It disagrees with the results 
of Wang et al. (1982) which range from 8®C day”l at the equator to only,13“C 
day“l at 90“S. The latter authors note their differences with Park and 
London, and ascribe It to different ozone distributions. However, at these 
altitudes and locations the ozone distributions are very similar. The 
possibility of excess ozone heating in the model cannot be discounted or 
confirmed given the uncertainties which exist for heating rates In the upper 
stratosphere. 

The temperatures calculated for January were Insensitive to the value of 
the drag used. In a separate experiment the drag was reduced by a factor of 4 
from that used In the standard run. Zonally averaged temperatures did not 
vary by more than 2®C at any location at any altitude after a one month 
Integration. 

JULY: Results for this month can be compared with the observations 

compiled by Oort and Rasmussen (1971) and Taljaard et al. (1969) for the 
troposphere, and by Groves (1971), Koshelkov (1977), and Labltzke (1980) for 
the stratosphere and lower mesosphere. The results shown by Hartmann (1976a) 
for July through September of 1973 in the Southern Hemisphere can be used to 
provide a check from the surface through the lower mesosphere. 

The simulation shown In Fig. 2 is In excellent agreement with the 
observations up to the 0. Imb level (65km) with the exception of the cold 
tropical temperatures from 50-60km, the same problem noted In January. 
Successfully reproduced features Include the very cold Southern Hemisphere polar 
temperatures In the lower stratosphere, the lower stratosphere temperature 
gradient reversal from the equator to 43**S, the temperature gradient reversal 
throughout the Northern Hemisphere lower stratosphere, and the magnitude of the 
polar temperatures In both hemispheres In the upper stratosphere and lower 



mesosphere. The region above 0. Imb is Influenced by the incorporated drag, as 
well as the model top, and thus the temperatures are not sufficiently cold in 
the Northern Hemisphere polar region. 

The lower stratospheric water vapor is not as excessive in polar regions as 
it was in January, due either to the extremely cold upper troposphere 
temperatures or simply the length of the run up to July as opposed to the 
following January. However, relative humidities are still greatly In excess of 
100^, and there would likely be an influence on the temperature if the moisture 
were allowed to condense. 

2) LATITUDINAL MEAN TEMPERATURE VERSUS MONTH 

Shown in Fig. 3 are the temperatures as a function of latitude and month 
for the 1.6mb layer (45km) representative of the upper stratosphere and the 34mb 
layer (23km) representative of the lower - to mid-stratosphere. These levels 
were picked to allow for comparison with satellite derived radiances centered on 
the 2mb and 30mb levels (e.g. , Labltzke, 1974, 1980, 1981; Labltzke and Barnett, 
1973; Ghazl, 1976; McGregor and Chapman, 1979). The 30mb temperatures are also 
obtained by radiosonde data and analyses are available (e.g., Labltzke and 
Goretzki, 1982; Naujokat, 1981), while at 45km additional observations come from 
rocketsondes (Groves, 1971; Koshelkov, 1977). 

In the model Northern Hemisphere, the upper stratosphere temperatures are 
coldest In December, warming somewhat In January, then cooling again in 
February, before the onset of the spring warming in March. The summertime 
temperatures show a peak around the solstice In June. In the lower 
stratosphere, high latitude northern temperatures remain cold throughout the 
winter. The relatively warmer January at 1.6mb results from warming pulses 
which occur during the month (Section 3f). These pulses do not occur with the 
same amplitude in the lower stratosphere, so it resembles the "cold” January 
observations. These features are in complete agreement with the generally 
occurring patterns noted in the references given above, although the 30mb polar 
temperatures are too cold by up to 10°C (Naujokat, 1981). 

In the tropical upper stratosphere a semiannual variation is apparent with 
coldest temperatures around the solstices. This is also In agreement with 
observations (Labltzke and Barnett, 1973) and has been correlated with 
temperature changes at high latitudes in the respective winter hemispheres 
(Fritz and Soules, 1970; van Loon et al. , 1972). As shown in Section 3f, this 
correlation occurs in the model as well. 

In the Southern Hemisphere temperatures are actually coldest in late fall 
in the upper stratosphere, and there is gradual warming during the winter. In 
the lower stratosphere the coldest temperatures do occur in midwinter. Both 
features are in accord with observations (Labltzke and Barnett, 1973; Labltzke, 
1974). 


Comparing the two hemispheres, in the upper stratosphere during midwinter 
(January versus July), the Southern Hemisphere temperatures are about 9**C 
colder than corresponding Northern Hemisphere temperatures. In the later part 
of the winter (February versus August) the situation has reversed, with the 
Southern Hemisphere now being 10®C warmer than the Northern Hemisphere poleward 
of 60® latitude. The warming pulses in the Northern Hemisphere in midwinter 
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Zonal mean temperature versus month for (a) 1.6mb level and 
(b) 34mb level. Dashed contours Indicate negative numbers. 
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are followed by intensified cooling in late winter. In contrast to the gradual 
warming scenario which prevails in the Southern Hemisphere. The magnitudes, 
locations and timing of these differences are all In agreement with satellite 
observations (Labltzke, 1974). 

The summer temperatures In the upper stratosphere at this level are 
warmer in the Southern Hemisphere In both model and reality, a fact which has 
been related to the shorter distance of the Earth from the Sun during Southern 
Hemisphere summer. In the model there Is a difference In solar radiation 
heating rates of 0.5-l“C day~^ between the two hemispheres during summer. The 
magnitude of the temperature difference Is In general agreement with 
observations (Labltzke, 1974; McGregor and Chapman, 1979). 

The Southern Hemisphere high latitude temperatures In the model's upper 
stratosphere are warmest In November; this type of spring "overshooting”, 
sometimes rapid, has been seen In the Northern Hemisphere (e.g. , Nordberg et 
al. , 1965; Johnson and Gelman, 1968). Theon et al. (1972) pointed out that 
the time of the temperature maximum in the Northern Hemisphere generally 
varies with latitude, occurring In spring at mid-latitudes, and In midsummer 
at the highest latitudes. This effect does occur In the model Northern 
Hemisphere and appears as a phase shift between the warm equinoctial 
conditions In the tropics and the warm solstice at the North Pole. In the 
model Southern Hemisphere conditions are warmest at the solstice from 
27®-67“S, and the November overshooting Is limited to the highest latitudes. 
The cause of this warming In the model remains to be Investigated, as does Its 
Interannual variability. 

In the lower stratosphere, fhe dominant hemispheric differences are at 
high latitudes, where the Southern Hemisphere Is colder In winter and warmer 
In summer, in agreement with what Is observed (Labltzke, 1974; McGregor and 
Chapman, 1979). The colder winters are thought to be due to the reduced 
planetary wave activity In the Southern Hemisphere, with the warmer summers 
again resulting at least partially from the shorter Earth-Sun distance. Both 
effects occur in the model (especially for standing waves - see Section 3d). 

3) LATITUDE-LONGITUDE TEMPERATURE MAPS 

The temperatures at various levels In the atmosphere are shown In Figs. 4 
and 5 for January and July. The top panels In each figure give the results 
for the surface, and for the mid-troposphere (mean temperature for the 
700-500mb layer, ~4.3km). The tropospheric results have been discussed in I, 
and will not be commented on further here except to note that they are In good 
agreement with observations. Results for the lower stratosphere (mean 
temperature for 100mb-30mb layers, ~20km) and mid-stratosphere (10mb-3.4mb, 
~35km) are given In the middle panels. The bottom panels show the results for 
the upper stratosphere (3.4mb-0.7mb, ~45km) and lower mesosphere 
(0. 7mb-0. 16mb, ~55km). 

JANUARY: The lower stratosphere temperature distribution features a warm 

region centered around the International Dateline In the Northern Hemisphere 
with a cold polar region centered around the Greenwich Meridian. Comparison 
with climatological lOOmb maps (Crutcher and Meserve, 1970) and 30mb maps 
(Labltzke, 1977, 1980; Labltzke and Goretzkl, 1982) Indicates that the 
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Fig. 4. Global distribution of temperature in January for (a) surface; 
(b) 700-500mb; (c) 100-30mb; (d) 10-3. 4mb; (e) 3.4-0.7mb; 

(f) 0.7-0. 16mb. Dashed lines indicate in-between contour values. 
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simulation Is In good agreement with observations. The polar temperatures are 
slightly too cold, as discussed earlier, but the difference Is small at 30mb 
compared to "cold" January temperatures (e.g., 1981 as shown by Labltzke and 
Goretzkl, 1982). 

In the middle stratosphere the axis of the warm and cold regions has 
sloped to the west, to a position In agreement with observations (e.g., the 
observations shown by Manabe and Mahlman, 1976; McGregor and Chapman, 1979; 
Austen et al. , 1977). The amplitude of the essentially wave 1 feature In the 
model upper stratosphere Is greater than usually observed, as will be 
discussed In Section 3d. In the lower mesosphere the warm region Is 
approaching the pole, as occurred In the observation of this altitude shown by 
Austen et al. (1977). This reversal of the latitudinal temperature gradient 
helps close off the stratospheric Jet (from the thermal wind relationship), as 
will be shown In section 3b. 

In the Southern Hemisphere during this month the temperatures are warmest 
near the pole throughout the stratosphere and lower mesosphere with a 
symmetric distribution (circular around the pole In polar stereographic 
coordinates; see Taljaard et al. , 1969 who show that this Is true even at 
lOOmb). 

JULY; The monthly average Southern Hemisphere winter has a much more 
symmetric temperature pattern than did the Northern Hemisphere, as Is expected 
from the reduced standing long wave energy (Section 3d). The warm ridge that 
becomes evident by the upper stratosphere Is located just to the west of the 
Dateline, In apparent agreement with the observations of McGregor and Chapman 
(1979) (but not with those shown by Labltzke, 1977 for June). The temperature 
gradient near the pole begins to reverse at heights of 45-50km, and Is 
completely reversed In the lower mesosphere, acting to close off the 
stratospheric jet. This appears to be an observed feature above Imb 
(Labltzke, 1980). 

In the Northern Hemisphere the warm polar temperatures and the concentric 
appearance of the Isotherms Is In accordance with observations (e.g., Crutcher 
and Meserve, 1970; Manabe and Mahlman, 1976; Labltzke, 1980). 

4) TEMPERATURE, LONGITUDE VERSUS ALTITUDE, 51 ®N 

Shown In Fig. 6 Is a longitude versus altitude plot of temperature at 
51®N for January and July. In January the warm regions start near the surface 
from what Is essentially a wave 2 feature at 10”W and ISO^W, slope to the west 
with Increasing altitude and become essentially a wave 1 feature In the 
stratosphere. The warm ridge Is located at about 150“E at 30mb and shifts 
westward to near 60“E at 1.6mb (43km), and to 60®W by 70km. These results are 
In good agreement with those shown by Labltzke (1980, 1981) for temperatures, 
or. In the upper stratosphere and mesosphere, radiances. 

The results for July show the absence of any long wave features above the 
tropopause, as expected due to the presence of east winds prohibiting wave 
propagation (Sections 3d, 3e). 



Pressure (mb) Pressure (mb) 


(a) 



quality 



West Longitude East 


Fig. 6 


Longitude-height profile of temperature at 51 "N In (a) January and 
(b) July. Dotted lines Indicate In-between contour values. 



b. Wind 


1) ZONAL MEAN WIND VERSUS ALTITUDE 

Shown in Fig. 7 are the results for January and July of the zonally 
averaged zonal wind versus pressure. 


JANUARY: The Northern Hemisphere shows the two prominent jets, 

tropospheric and stratospheric. The tropospheric jet is slightly weaker than 
observed (in I it was shown that it improves in intensity when the horizontal 
resolution is Improved). The stratospheric jet appears realistic (e.g. , 
Newell, 1968) although once again large variations are encountered in 
observations due to the occurrence of stratospheric warming events with light 
west winds. It will be shown below that the value of the jet is actually 
slightly larger in December and February, since January underwent warming 
pulses in the stratosphere. When the drag in the top three levels was reduced 
by a factor of four the stratospheric jet Increased by 5-10ms“^. 

The vertical shear to the north and above the tropospheric jet is too 
weak in the model, a feature common to models with topography (Tenenbaum, 
1982). This is consistent with the cold polar lower stratosphere discussed 
above. In his analysis Tenenbaum concluded that this deficiency affected wave 
3, allowing it to propagate too easily into the stratosphere, and thus 
reducing the tropospheric wave 3 energy below its observed value. As will be 
shown in Section 3d, the effect may occur in this model as well. 


Observations of winds in the extratroplcal mesosphere often give values 
of 70 ms“^ near 75km (Newell, 1968; Groves, 1969), although there are many 
observations of significantly smaller velocities (e.g., Gregory et al. , 1981). 
The model value at that altitude in January is 45 ms"^ at upper mid-latitudes. 
While the observations will be affected by whether or not warmings have 
occurred, it is possible that the drag incorporated in the model, which is 
responsible for the decrease of wind speed above the stratospheric jet, might 
be too large. 


The low latitude easterlies in January cross over into the Northern 
Hemisphere. In Section 3e it will be shown that this is related to angular 
momentum divergence caused by the long waves in the winter hemisphere, as 
postulated (Hlrota, 1980). South of the equator the low latitude east winds 
appear reallsTlc up to about 45km, but then become excessive. The geostrophlc 
wind shear is very sensitive to the temperature gradient near the equator. 

For example, a 1“C temperature gradient between 12® and 4® latitude at 1.6mb 
will produce a wind shear of about 11 ms“^ up to 0.7mb; thus if the gradient 
Increases to 3®C, the wind shear will be 33 ms”^ (which is similar to what 
occurs at low latitudes in the model). Given the potential problems with the 
radiation at 50km and above, the excessive temperature gradient right near the 
equator is not surprising. A similar problem was encountered by Crane et al. 
(1980) in their two dimensional model. 


The direction of the zonal winds in the tropical lower stratosphere 
depends upon the phase of the quasibiennial oscillation. Given the length of 
this integration it is impossible to tell whether the model will experience 
such an oscillation. During the model run the winds in this region were 



Pressure (mb) Pressure (mb) 




7, Latitude-height profile of zonal mean zonal wind for (a) January 
and (b) July* West winds are shown by solid contours , east winds 
by dashed contours. 





generally easterly although a region of west winds at the equator descended 
from 73mb to ISOmb during the course of the run. 

The Southern Hemisphere simulation, with tropospheric westerlies and 
stratospheric east winds appears realistic, although the subtropical 
tropospheric Jet stream may be too weak (Labltzke, 1980, Fig* 5). 

JULY: The July simulation correctly reproduces the summer stratospheric 

east winds and the winter stratospheric Jet. Note the stronger magnitude of 
the Southern Hemisphere Jet compared to that of Northern Hemisphere winter in 
the stratosphere; the values are in agreement with observations (McGregor and 
Chapman, 1979). Once again the tropospheric Jet streams may be slightly too 
weak, and the east winds to the north of the equator are excessive above 50km. 
Again note the region of decreased east winds, this time manifesting as actual 
west winds, at 1.6mb, now Just to the south of the equator. The mesospheric 
. east winds in the Northern Hemisphere, which are also affected by the magnitude 
of the drag, are not too strong compared to observations (Groves, 1969). 

MERIDIONAL VTIND: Fig. 8 shows the meridional wind for January and July. 

It is difficult to compare these results with observations in the 
stratosphere; meridional winds in summer are weak as the circulation is 
symmetric around the pole, while in winter the wind field is dominated by 
large scale waves, and it would be necessary to properly incorporate 
observations from all longitudes to get a meaningful average. In addition, 
tidal variations near the stratopause are probably larger than the prevailing 
wind, so rocket observations, often taken near noon, introduce a spurious 
tidal component. The model tended to produce equatorward flow in the winter 
hemisphere at latitudes poleward of about 50“ throughout the stratosphere. 

This appears to be in general agreement with the analysis of Crane et al. 

(1980) for December, 1976. 

In the mesosphere the circulation from the summer pole to the winter pole 
is in agreement with observations (Groves, 1969) and model results (Holton and 
Wehrbeln, 1980), although it is perhaps at a slightly lower altitude. This 
circulation is an artifact of the drag included in the top three levels, as it 
will act on both east and west winds to produce summer to winter pole flow. 

The magnitude of the flow is excessive for this altitude, reaching 10 ms“^ in 
the upper layer, while the (admittedly uncertain) observations give only a few 
ms“^. At 90km, however, values of 10-15 ms have been found (Nastrom et al. , 
1982). In the lower mesosphere below the layers with drag the circulation 
produces south winds from 20“S to 40“N, and north winds elsewhere - this is in 
agreement with observations (Groves, 1969). 

INERTIAL INSTABILITY; The meridional wind in the tropical stratosphere 
shows what appears to be inertial instability. The evidence is as follows: 
the meridional wind oscillations in the vertical appear clearly in both the 
months shown in Fig. 8 in the hemisphere in which the criterion for inertial 
instability (f-dU/dy less than zero - Northern Hemisphere; greater than zero 
in the Southern Hemisphere) is exceeded - Northern Hemisphere in January, 
Southern Hemisphere in July. As noted above, the east winds are found across 
the equator in the winter hemisphere, in both observations and the model, due 
to the influence of planetary wave momentum divergence at low latitudes of the 
winter hemisphere. The inertial instability extends only up to the latitude at 
which the above criterion is actually met; dU/dy is of order 2-4 x 10~^ s~^ a 
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Fig. 8. Latitude-height profile of zonal mean meridional wind for 
(a) January and (b) July. South winds are shown by solid 
contours, north winds by dashed contours. Off-contour numbers on 
the graph Indicate values In specific locations. 





1.6mb In the model, and It drops below the value of the Coriolis parameter at 
about 20° latitude. During the equinoxes when cross equatorial shear is 
largely absent, these oscillations are either completely gone or substantially 
reduced. The oscillations are apparent in all the fields of interest, 
including the vertical velocity and thus the stream function, and eddy 
transports (Section 3e). 

This oscillation appears to be similar to that obtained by Hunt (1981) In 
a general circulation model experiment. According to the analysis of Dunkerton 
(1981), cross-equatorial shears on the order of 1.5 day”^ should have vertical 
wavelengths on the order of 1km, and thus not be resolvable in models with 
coarser vertical resolution* The shears in this model were on the order of 2.3 
day“^, which is consistent with the greater vertical resolution used (5km), 
although shears about three times as large would be necessary to be 
quantitatively consistent with Dunkerton* s analysis. The coarse resolution, 
both vertical and horizontal, may well exaggerate the magnitude of this effect 
in the model. 


The vertical velocity perturbations induced by this instability, on the 
monthly average, are on the order of 5 x 10~^ mb s“^. If this were to go 
entirely into local heating it could produce heating rates of 2“C day“^ at 
1.6mb, which is certainly significant. The temperature profile does not show 
obvious perturbations in this region (Fig. 2), but what the effect of any 
heating is on the latitudinal temperature gradient and thus the zonal wind 
shear remains unknown. 


The zonal wind field also shows perturbations (Fig. 7).‘ These are located 
at or near the equator, in the region which has been characterized as 
associated with inertial instability. The zonal wind perturbation even 
switches hemispheres with the Inertial instability, and is thus likely 
associated with it. However, the magnitude of the zonal wind oscillation is 
very large, with peak to peak values of 50 ms“^, centered on the equator at 
1.6mb. This is an order of magnitude larger than the meridional wind 
oscillation, and much larger than the oscillation shown by Hunt (1981) in the 
zonal wind. It might be associated with heating gradients Induced by the 
vertical velocity perturbations. It is associated with local regions of 
angular momentum convergence and divergence and might be related to the effect 
of the instability on propagating Kelvin waves. Kelvin waves are discussed 
further in Section 3d. 


2) LATITUDINAL MEAN ZONAL WIND VERSUS MONTH 


Fig. 9 shows the variation of the zonal wind at 1.6mb and 34mb as a 
function of latitude and month. In the Northern Hemisphere, the upper 
stratosphere picture indicates that west winds dominate during the winter, 
peaking in February and being slightly weaker in January than the other winter 
months. The model's warming pulses in January, similar to those often observed 
in the real stratosphere during particular winter's, were responsible for the 
weakening in January. In the lower stratosphere there is a more progressive 
Increase of west winds with time, again peaking in February. Summer easterlies 
appear at both levels, stronger in the upper stratosphere* This simulation is 
in agreement with observations at Northern Hemisphere extratropical latitudes 
(e.g. , Richards, 1967; Groves, 1971) except for the winds in the lower 
stratosphere winter being excessive by about 10 ms“^, as discussed earlier. 
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(The east winds at both poles at both levels are a numerical artifact which 
results from using the Coriolis force and metric term at the pole). 

In the Southern Hemisphere the seasonal variation Is similar, except that 
the winter west winds are stronger than was the case during Northern Hemisphere 
winter. This Is an observed feature (McGregor and Chapman, 1979 for the 
geostrophic wind field) and results partly from the lack of major warmings to 
Interrupt the west wind field during winter. The lower stratosphere west winds 
are not excessive, as the cold polar stratosphere was properly simulated. 

SEMIANNUAL WIND: The tropical upper stratosphere Is dominated by a 

semiannual wind oscillation (e.g. , Belmont et. al., 1974). Due to the zonal 
wind perturbation associated (In some sense) with the Inertial oscillation, the 
1.6mb level shown in Fig. 9 is not representative of the tropical upper 
stratosphere for this parameter. The oscillation which is present during the 
solstices produces light east or west winds at this level. In contrast to the 
observed flow (and flow at other levels) which shows maximum east winds during 
these seasons. If one averages the wind at the equator between 3.4mb and 
0.17mb (38-60km) one obtains the variation shown in Fig. 10. East winds now do 
peak during the solstices, but the west wind phase is weaker than observed 
(Reed, 1966) (and in fact the April/May winds show only a diminution of east 
wind intensity). If the west wind phase is due to the momentum transport by 
propagating Kelvin waves (Hirota, 1980) this Implies that either the Kelvin 
wave generation is too weak, or the waves are not reaching the upper 
stratosphere. This is discussed further In Section 3d. 

3) ZONAL WIND, LONGITUDE VERSUS ALTITUDE, 51 *N 

Fig. 11 shows the variation of the zonal wind versus longitude at 51 “N in 
January and July. The winter picture shows dramatically the dominant wave 1 
feature In the stratosphere, with Increasing wave number In the troposphere. 

At this latitude the strongest tropospheric winds occur in the eastern Aslan 
region (Lau, 1978) as simulated by the model, and slope westward with height 
to become associated with the west wind maximum at upper stratospheric levels. 
As will be shown, the phase of wave 1 appears proper, although the amplitude 
In the upper stratosphere Is probably excessive. 

In July, the stratospheric westerlies are replaced by east winds, there 
Is little wavff" propagation, and the stratospheric wind field shows little 
longitudinal variation. 

c. Geopotential height 

Figs. 12 and 13 show the geopotential height fields for various levels In 
the troposphere and stratosphere for January and July, as well as the sea 
level pressure. 

JANUARY: The model produces realistic subpolar lows In both 

hemispheres, although the low centered at 30“E near Antarctica may be too weak 
(Taljaard et al. , 1969). There is also a tendency for the subtropical highs 
to be displaced poleward of their observed positions. As shown in I, both of 
these deficiencies are removed by using greater horizontal resolution; in 
addition, as also shown In I there Is much year to year variability In sea 
level pres-sure In the model, so It Is improper to fully evaluate this 
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Pig. 10. Zonal Bean equatorial wind averaged from 3.4mb-0.I7mb versus month. 
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Longitude-height profile of zonal wind at 51 ®N in (a) January and 
(b) July. Dashed contours indicate east winds, dotted contours 
indicate in-between contour values. 
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Fig* 12. Global distribution In January of (a) sea level pressure, and 
geopotential height fields at (b) 500mb; (c) lOOmb; (d) lOnb; 

(e) 3.4nb; (f) 0.7nb. Sea level pressure Is In mb with contours 

every lOmb; all others In hectometers (100m). Dashed lines represent 
In-between contour values. 
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Fig. 13. Same as Fig. 12 for July. 











parameter on the basis of a one year simulation. At 500mb the Northern 
Hemisphere simulation looks reasonable both in phase and amplitude (compared 
with Crutcher and Meserve, 1970) except that the short wave feature centered 
in the model at 60®E in mid-latitudes may be slightly too far east. The 
comment on Interannual variability applies at this level as well, and this 
particular feature has large observed variability (Lau and Oort, 1981). In 
the Southern Hemisphere the lack of a strong standing wave pattern, along with 
uncertain observations, makes comparison difficult, although the trough/rldge 
feature shown over Australia appears accurate (Taljaard et al., 1969). The 
tropospheric simulation is important as it establishes the position and 
amplitude of the standing waves which propagate into the stratosphere. 

By lOOmb in the Northern Hemisphere the observations show troughs 
centered over North America, Europe and Japan, with ridges in the Atlantic and 
Pacific. The model reproduces these features although the European trough is 
again too far east and probably too weak. As we will show in the next section 
this represents an underestimate of wave 3 energy in the troposphere. In the 
Southern Hemisphere the polar low is correctly positioned in the eastern 
hemisphere. 

At lOmb the Aleutian High is properly positioned and at approximately the 
correct height (e.g. , c.f. Lin, 1982), while the polar low is about 1km too 
low, due to the cold lower stratosphere temperatures. Climatological 
observations for the Southern Hemisphere for this month are difficult to 
obtain; the symmetric circulation around the pole is apparently correct 
(Cunnold et al. , 1975), due to the inability of tropospheric waves to 
propagate vertically in an east wind regime, while the computed height of the 
polar anticyclone is slightly greater than that computed by Manabe and Mahlman 
(1976). 

At higher levels, the Aleutian High amplifies and shifts slightly to the 
west. Perusal of the weekly synoptic analyses at these levels (Staff, Upper 
Air Branch, 1967-1971) shows that this is a realistic feature; the comparison 
between the mode wave 1 amplitude and that of the observations will be given 
in section 3d. In the Southern Hemisphere the circulation remains 
symmetric into the lower mesosphere. 

JULY: Fig. 13 shows the geopotential height surfaces for this month. 

The subpolar lows in the Southern Hemisphere are of approximately the right 
magnitude, although their position shows some variance with climatology 
(Taljaard et al., 1969). The subtropical highs in both hemispheres are too 
far poleward. At 500mb the Northern Hemisphere has shorter wave features, as 
is expected for the season, but the subtropical highs are not well formed. In 
the Southern Hemisphere the weak troughs adjacent to the continents appear to 
be consistent with observations. By lOOmb the model shows high pressure right 
near the poles: as mentioned in Section 3b this is affected by the numerical 

formulation at the poles. Northern Hemisphere observations (Crutcher and 
Meserve, 1970) show weak low pressure near the pole. The major deficiency is 
again the lack of well formed subtropical highs. In the Southern Hemisphere 
the simulation seems realistic, with the polar low displaced towards the 
Austrailian sector. 

At lOmb the Northern Hemisphere circulation is properly symmetric about 
the pole (Staff, Upper Air Branch, 1967, 1969) and it remains so into the 



lower mesosphere* ‘ In the Southern Hemisphere the winter circulation Is also 
fairly symmetric, around the polar low. Indicative of the dominant wave number 
zero pattern In the standing geopotential field (e.g. , Labltzke, 1977). The 
long wave amplitudes In both hemispheres will be examined In Section 3d. 

d. Eddy energy 

1) EDDY KINETIC ENERGY 

Fig. 14 shows the eddy kinetic energy as a function of latitude and 
pressure for January and July. The distribution In the Northern Hemisphere 
troposphere Is realistic In both seasons except that In January the eddy 
energy Is somewhat excessive at low levels, while In July It Is slightly 
deficient In the upper troposphere. In the Southern Hemlsphere^values appear 
realistic In both seasons. This Is essentially the same conclusion given In 
I, to which reference should be made for details of the comparison with 
observations. 

In the equatorial tropopause region In July a local excess of eddy energy 
appears. This Is essentially a wave 1 standing feature, which derives Its 
energy from potential energy, generated by latent heat release. It Is not at 
all evident In June, and la greatly diminished in August; nor does It exist In 

the five year average shown In I. This feature may be caused by an Inertial 

Instability as wind shears of the right sign do exist, but Indications are not 
as clear as for the oscillation near the stratopause* 

The eddy energy (per unit mass) increases with altitude In the winter 
hemisphere, much more so in the Northern than Southern Hemisphere. The summer 
hemispheres show little eddy energy in the stratosphere. These are all 
reasonably realistic features; It Is difficult to compare with observations 
which are not usually given In this form In the stratosphere. A better 
comparison will be possible In the sections below. 

Evident In both months is a local eddy energy maximum In the stratopause 

region near the equator, possibly associated with the Inertial Instability. 

2) ENERGY CONVERSIONS AND BUDGET 

How Is the energy generated In the stratosphere? Figs. 15 and 16 show 
the energy diagrams for the troposphere and several regions In the 
stratosphere and lower mesosphere during January and July. These differ from 
the usual presentation of such diagrams In two ways. First, the wide arrows 
entering or leaving the zonal and eddy kinetic energy boxes represent the 
effect of the full pressure gradient force term. Including both the 
convergence of geopotentlal flux and the barocllnic conversion from potential 
to kinetic energy; a decomposition of the two terms with better vertical 
resolution will be shown below. Second, due to the nature of the diagnostic 
routine, the transformation from zonal available potential energy to eddy 
available potential energy could not be calculated. 

In both months the winter hemisphere has more zonal and eddy kinetic 
energy up through 0.48mb than the summer hemisphere. The Northern Hemisphere 
has greater eddy energy, both kinetic and potential in winter than does the 
Southern Hemisphere, while the reverse is true for the zonal potential and 
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Fig. 15. Energy budget for January In varloua regions of the atmosphere. 

Schematic shown at top, with zonal available potential energy (AZ), 
eddy available potential energy (AE), zonal kinetic energy (KZ) and 
eddy kinetic energy (KE), In units of lO^^J. Transformations shown 
include KZ to KE, radiative generation or dissipation of AZ and AE, 
dissipation of KZ and KE, and effect of pressure gradient force 
(broad arrow) on KZ and KE (see text) In units of .^g 

various effects on A and K In the troposphere are not shown - see I 
for more details. Northern Hemisphere on left. 





July 


Rad — PISS 21 Layer 8°xl0‘ 


PG 

Force 


Rad 


— I 


DISS 


PG 

Force 


N.H. 


Lower Mesosphere 
0.48-0.03 MB 


OH 


0.4 




6.9 

0.3 

iL 


26 h6.5 


1.8 


S.H. 

0.3<f^ y^fliT^II 

8.7 T 
3.3 

i 


2.2 b 1.9 


0.34^0.7 ^7.9b4.4 


3.3^30 I ^|Tl9l 
3.3 

1.2 


Mid and Upper Stratosphere 

10-0.48 MB 42, 


103 


0.2< l30l 




302 

T 

36 


5<rW\ 

— 39^^ 




Lower Stratosphere 
100-10 MB 


r 

0 

i 


3C^ 

I- Hi022| '^|9I6 


30 


9« 43701 ^H05 

0^ 


I2<^ 253 


30 


360 


[iiiel ^^JioT 


Troposphere 
984-100 MB 


21313 


360 


1 4973 


246 


30 


1453 


4512 


2567 

“T" 

114 

_L 


-<5213019 
264 


Fig. 16. Sane as Fig. 15 for July. 




kinetic energy. The Southern Hemisphere Is dominated by a very cold polar 
low, while the Northern Hemisphere has large standing waves In winter. Both 
potential and kinetic energy decrease with altitude along with the mass of the 
atmosphere. 

The magnitudes of the zonal and eddy potential and kinetic energies are 
in rough agreement with observations for specific months (generally 
±50% - Dopplick 1971, Hartmann, 1976b) although comparisons are difficult since 
neither study Included tropical areas, which, area weighted, occupy one-third 
of the globe. 

The tropospheric energy cycle Indicates that eddy kinetic energy Is being 
generated by the pressure gradient term. The magnitudes of the transformation 
agree well with observations (e.g., Oort, 1964) as noted In I. In both winter 
stratospheres the energy cycles also show that the pressure gradient term Is 
generating eddy kinetic energy, which In these regions Is completely 
transferred Into zonal kinetic energy, and then completely transferred out. 

Into zonal potential energy or as a vertical flux. This Is Indicative of an 
energy cycle driven by waves propagating through the region, rather than 
In-situ processes. (See Plumb, 1983 for an alternate version of the energy 
cycle, which clarifies this point). Radiative fluxes destroy eddy available 
potential energy and generate zonal available potential energy. The energy 
transformations are again In rough agreement with observations (Dopplick, 

1971; Hartmann, 1976b). 

The summer stratospheres show, in general, smaller energies (except for 
the zonal available potential energy In January), and negligible energy 
transformations or fluxes. 

The mesospheric energy cycle Is strongly Influenced by the large kinetic 
energy dissipation which occurs in the model at those levels. To the extent 
that this dissipation also arises In the real atmosphere. It must be replaced 
by processes such as generation from the pressure gradient force, as in the 
model. 


The decomposition of the pressure gradient force Into the eddy generation 

term (-m'o*) and the eddy convergence term (-3/3p (u)'^') - v' • V^') is 
shown In Fig. 17 for both months and hemispheres. In the troposphere energy 
Is generated barocllnlcally between about 300mb and 900mb, accompanied by 
transport out of this region (thus a net divergence), both upward and 
downward. The upward flux propagates Into the stratosphere (especially In the 
winter hemisphere extratropics) and Its convergence provides the kinetic 
energy source indicated in Figs. 15 and 16. The convergence term thus 
dominates the eddy kinetic energy generation term In the model winter 
stratosphere, in agreement with observations (Newell and Richards, 1969; 
Miller, 1970; Dopplick, 1971; Hartmann, 1976b). As the energy Is, to a good 
extent, propagating vertically straight through the region (to the extent that 
the noninteraction theorem conditions are obeyed) the energy simply cycles 
Into eddy kinetic energy, then to zonal kinetic energy, and then out, as shown 
In Figs. 15 and 16. 

The eddy conversion term Is generally opposite in character to the 
convergence term, not only on the hemispheric average but also at each 





Fig. 17 


Eddy generation (-w'o' ) and convergence (-3/3p Tu^VT ” v' • V4>') 
per unit o as a function of altitude In the Northern and Southern 
Hemispheres for (a) January and (b) July. 


latitude and altitude. The latitude-height pictures of conversion and 
convergence from Che model (not shown) look extremely similar to those shown 
by Manabe and Mahlman (1976) from their GCM. As shown In Fig. 17 the 
generation term Is slightly negative throughout the stratosphere; for steady 
waves In an adiabatic atmosphere both the vertical energy propagation and the 
barocllnlc conversion should be positive when the northward transport of 
sensible heat by the eddies Is downgradlent (thus In a region where west winds 
are Increasing with altitude) (Hartmann, 1976b). As we will show In Section 
3e, the model eddies do transport heat northward In the stratosphere where 
there Is positive shear, and yet the barocllnlc conversion Is slightly 
negative. This Is also the result found In both models and some observations 
(Manabe and Mahlman, 1976, Fig. 9.4; Kasahara and Sasamorl, 1974; Manabe and 
Hunt, 1968; Oort, 1964) although not by Doppllck (1971) or Hartmann (1976b). 

A possible reason for the discrepancy Is that the conditions are not 
adiabatic, and the phase relation between the temperatures and heights are not 
necessarily Chose of a steady or freely propagating wave. For example. If one 
adds a dissipation term to equation (6) In Hartmann (1976b), and uses the 
dissipation time constant associated with Newtonian cooling Indicated by 
Dickinson (1973), one finds that at 16mb the ratio between dissipation and the 
model's northward transport of sensible heat by eddies Is equal to 2.6 at 51 °N 
In January. Thus the sensible heat transport Is not balanced primarily by 
vertical motion, as Is the case for a free wave propagating In an adiabatic 
atmosphere, but Is dissipated effectively by radiative damping. 

The model stratosphere, and presumably the real stratosphere. Is made up 
of waves propagating with greater or less difficulty, depending upon 
dissipation, reflection or other phenomena (such as critical layers). For a 
vertically propagating wave, for which the vertical velocity and geopotentlal 
fields are In phase. If the temperature wave lags the height wave by more than 
1/4 wavelength It will result In a negative generation of eddy energy (as the 
warm air will be closer to the descending motion In the trough). However, the 
wave will continue to transport heat northward until the temperature wave lags 
the height wave by 1/2 wavelength (as, until that lag, the positive meridional 
wind will still occur In the warm sector). If one compares Fig. 4 and Fig. 12 
one can see that the standing temperature wave lags the standing height wave 
by approximately 1/4 wavelength (90® for wave 1), Indicative of little 
barocllnlc energy generation but maximum northward sensible heat transport for 
vertically propagating waves. The monthly average negative energy generation 
Indicates tha^ the average phase lag for all waves Is slightly greater than 
1/4 wavelength In the model stratosphere. An example of the day to day 
variation In this phase lag will be given In Section 3f. 

At levels above 2mb, the shear reverses, the northward heat transport Is 
now associated with a substantial negative generation term, and the eddy 
convergence Increases as wave dissipation becomes more effective. Despite the 
slightly negative generation term In the model stratosphere, compared to the 
slightly positive value found by Hartmann (1976b), the results In Fig. 17 look 
very similar to Hartmann's analysis of observations. 

3) SPECTRA AND LONG WAVE AMPLITUDES AND PHASES 

The spectra of the eddy kinetic energy averaged over various layers In 
the model are shown for January and July In Figs. 18 and 19. In both months 
the ratio of the short wave to long wave energy decreases as expected with 



Energy (lO'^J) 



Fig. 18. Eddy kinetic energy (lO^^J) as a function of wavenumber in January 
for various regions of the atmosphere. Northern Hemisphere on left. 
Note the change of scale with elevation. 
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Increasing elevation until the model has primarily wave 1 with some wave 2 
energy at the highest levels during winter. In summer there is little wave 
energy above the lower stratosphere. In the Northern Hemisphere winter wave 1 
energy appears to pass through the lower stratosphere relatively unaffected, 
and accumulate at higher levels, as Indicated for all waves in Fig. 17. In 
contrast, wave 2 energy (or activity) propagates in a more limited fashion 
between the same regions (it is reflected more easily). During July both wave 
1 and wave 2 show little energy throughout the stratosphere. In the 
mesosphere, wave energy is affected by the parameterized drag. 

The magnitudes and general shape of the spectrum for January in the 
Northern Hemisphere troposphere compare extremely well with the observed 
spectrum shown by Chen (1982) for the winter of 1976-77; the only deficiency 
is that the energy dip shown in the figure for wave numbers 3 and 4 was not 
observed during that winter. The January simulation which was part of the 
three month spln-up also had deficient energy in these wave numbers, but the 
alternate January run %d.th corrected water vapor had 30% more wave number 3 
energy. 

The somewhat smaller spectral elope seen through the medium scale wave 
numbers during Southern Hemisphere %d.nter in the troposphere is hard to 
verify, although the relative Importance of the medium scale waves during 
Southern Hemisphere summer has been documented (Randal and Stanford, 1983), 
and also appears In the model. Both van Loon and Jenne (1972) and Trenberth 
(1980) show that standing waves 1-3 account for most of the standing wave 

energy at SOOmb - this Is true in the model as well, as waves 1-3 in the 

Southern Hemisphere troposphere during July represented 60% of the standing 
wave energy. However, in the model the standing waves accounted for only 23% 

of the total eddy energy, and the spectrum shown in Fig. 19 fails to reveal 

the long wave dominance when transient energy is also considered. In the 
Northern Hemisphere troposphere during winter, standing wave energy accounted 
for 40% of the total eddy energy. 

More definitive comparisons can be made by looking at the stationary long 
wave geopotential amplitudes and phases - these are shown In Figs. 20 and 21 
for January and July for wave numbers 1-3. 

JANUARY; The Northern Hemisphere wave number 1 amplitude reaches very 
large values In the upper stratosphere then starts decreasing; the decrease 
begins below the region of the parameterized drag (which occurs above 0. Imb 
and is not shown In these figures), and appears to be caused by transfer of 
energy into zonal kinetic energy and radiative damping. The phase of wave 1 
slopes westward with height, fairly rapidly in the troposphere and lower 
stratosphere, and then more slowly above lOmb. The rapid phase progression Is 
Indicative of a relatively freely propagating wave, as noted in the discussion 
of the spectral diagnostic (Fig. 18), while the slower variation Implies less 
wave transmission. Wave 2 reaches a more modest amplitude maximized at a 
lower elevation (lOmb), and Its phase variation Implies little propagation 
above that height, as did the spectral analysis. Wave 3 grows to a smaller 
amplitude maximum In the upper stratosphere, and the spectral analysis shows 
that this Is associated with substantially reduced energy; Its phase variation 
indicates little propagation above the lower stratosphere. The relative 
amplitudes and phase variations among these three waves are In accordance with 
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Fig. 20. Latitude-height profile of the amplitude (m) and phase (location of 
ridge, measured west of Greenwich) in January for wave numbers 1, 2, 
and 3 (from top to bottom). For the phase of wave number 1 and all 
wave amplitudes, dashed lines indicate in-between contour values; for 
the phase of waves 2 and 3 they indicate degrees east longitude. 
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expectation given the increased susceptibility to west wind reflection 
encountered by waves with increasing wavenumber (Charney and Drazin, 1961). 

There are numerous observational data analyses with which to compare the 
model results. As stratospheric dynamics rely heavily on long wave generation 
in the troposphere we will discuss the tropospheric simulations in some 
detail. The Northern Hemisphere January wave 1 amplitude of 100m in the 
mid-troposphere at 50“N, rising to 130m in the upper troposphere is in 
agreement with observations (Muench, 1965; van Loon et al. , 1973; Defant et 
al. , 1979; Reiter and Westholff, 1981; Geller et al., 1983). The phase 
variation, such that the wave starts near 50°E at 50°N near the surface, tilts 
westward so that the ridge is near the Greenwich meridian by 500mb, and near 
lOO^W by lOOmb, is also in agreement with these observations. At SO^N wave 
number 2 amplitude is about 40% deficient in the mld-tropospherg, although its 
value of 180m at lOOmb is in agreement with observations. The phase variation 
at this latitude, starting at 60®E near the surface, tilting to near 30“E by 
500mb and near 10®E by lOOmb is slightly greater (by about 10“ ) than that 
indicated by most of the observations, consistent with the more rapid growth 
experienced by the wave in the model. Wave number 3 is about 50% deficient in 
amplitude throughout the troposphere, as indicated in the spectral analysis, 
although its slope, starting near 0“ at the surface, and sloping just slightly 
to reach 30“W by lOOmb is correct. Thus the model provides proper amplitude 
wave 1 and wave 2 forcing to the lower stratosphere, which is highly Important 
as they are the predominant factors in driving stratospheric dynamics. The 
wave 3 amplitude, however, is deficient by a factor of 2; this should have 
little effect on the mid and upper stratosphere due to the difficulty this 
wave experiences in propagating into the stratosphere. (However, as shown 
below wave 3 amplitudes at higher levels in the stratosphere appear realistic. 
This lends some credence to the conclusion of Tenenbaum (1982) that wave 3 in 
the troposphere is underestimated in models because too much energy propagates 
into the statosphere due to the lack of sufficient negative shear above the 
jet stream level - and thus Incorrect refractive index. Uncertainty in the 
effect of photochemical dissipation, as is discussed below, prevents any firm 
conclusion). Wave 3 amplitudes are 33% larger in the mid-troposphere in the 
alternate January with corrected stratospheric water vapor, which again raises 
the question of model, and natural, variability. 

At lOmb the model has the following amplitude and phase (at 65“N, the 
location of m^lmum amplitude): wave 1 - 1000m, 170“E; wave 2 - 400m, 175“E; 

wave 3 (at 51 “N) -90m, 100“E. Comparison with observations (van Loon et al. , 
1973; Barnett, 1980; Quiroz, 1980; Geller et al. , 1983) shows that the 
amplitudes are all in agreement with observations, although both wave 1 and 
wave 2 amplitudes reach those values only in given years (1965, 1977, 

1978 - van Loon et al., 1973; Quiroz, 1980). The phases are all in good 
agreement. By Imb, observations (Barnett, 1980; Smith, 1983; Geller et al., 
1983) show amplitudes for wave 1 equal to, and for wave 2 less than those at 
lOmb, while the model shows wave 1 amplitude having doubled, and wave 2 
remaining relatively constant. Even considering the uncertainties in the 
observations (both Green, 1972 and Hirota and Barnett, 1977 report large 
amplitude standing waves in the mesosphere) it is likely the model 
overestimates the long wave amplitudes, for wave 1 and possibly 2, by 50% or 
more at Imb. Wave 1 energy at this altitude in the model is lost primarily 
through interaction with the mean flow (KE to KZ) and secondarily through 
radiative damping. However, the radiative damping provides a violation of the 



noninteraction theorem, and thus may Indirectly allow a portion of the energy 
transfer to the mean flow to occur. As noted in Sections 2b and 2d, the model 
does not Incorporate photochemical damping associated with the temperature 
dependence of the ozone reactions. Ghazi et al. (1979) have estimated that 
Inclusion of this damping Increases the total damping In the upper 
stratosphere by about a factor of two. Whether this would be sufficient to 
reduce the wave amplitudes to values more In line with published observations 
is uncertain. Wave 3, which was underestimated at lOOmb, reaches an amplitude 
of about 130m between 1 and 2mb, In agreement with observations (Geller et 
al., 1983), and thus appears to be of realistic amplitude throughout most of 
the stratosphere. Of course Its amplitude might also be reduced, at least In 
the upper stratosphere, were photochemical damping to be Included. Diffusion 
would also reduce the amplitudes of all the waves. The phase variation of all 
three waves between lOmb and Imb are In agreement with observations, wave 1 
shifting slowly westward to 160"E by Imb at 65°N, while waves 2 and 3 show 
little variation with altitude. As the phases of the waves are correct 
throughout the stratosphere, the waves appear to be propagating in a realistic 
fashion. 

In the Southern Hemisphere mid-troposphere during January, the model's 
values for wave 1 amplitude of 50-70m are slightly less than the observed 80m 
(van Loon and Jenne, 1972; Trenberth, 1980) at 50-60“S while the phase, near 
180° longitude. Is some 30° west of the observed phase. For wave 2 the 
amplitude of 20-30m In mid-latitudes appears realistic (Trenberth, 1980), 
while the phase is about 20° too far west; the higher amplitude of 50m at 67°S 
has not been reported. The wave 3 amplitude of over 60m Is more than twice 
that observed, while the phase of about 30°E from 50-60°S Is appropriate. In 
both the model and observations (van Loon and Jenne, 1972) wave 1 shifts 
westward with height at high latitudes up to about 30mb, and then the phase 
becomes stationary, implying that upward propagation effectively ends at 
that level (due to the appearance of the summer east winds as In Fig. 7a). The 
amplitude maximum of 140m at lOOmb is greater than that reported (from one year 
of observations In van Loon and Jenne, 1972). Wave 3 shows little phase 
variation with altitude and thus little propagation In either the model or 
observations. If the observations are accurate and representative, the model 
seems to overestimate the standing wave amplitude of wave 3 in the Southern 
Hemisphere troposphere, and wave 1 In the lower stratosphere, while maintaining 
an approximately accurate phase. The Southern Hemisphere summer stratosphere 
has very small wave amplitudes, as Is observed. 

JULY: The results for the different wavenumbers are shown In Fig. 21. 

In the Southern Hemisphere winter the tropospheric amplitudes for waves 1-3 all 
reach 100m, Increasing to 800m for wave 1 by 0. 7mb, to 250m for wave 2 at 
3.4mb, and with a slow but general decrease with altitude for wave 3. In the 
Northern Hemisphere mid-stratosphere, more than half of the eddy energy 
associated with wave 1 is in the standing waves, and slightly less than half 
for wave 2. In the Southern Hemisphere, only about 1/4 of each Is in the 
standing wave component - this agrees with observations, especially for wave 2 
(Hartmann, 1976a) (and as noted earlier,, is true for the Southern Hemisphere 
troposphere as well). This makes comparison with observations difficult above 
about 30mb for the standing waves, since analysis of continuous time series are 
necessary to remove the traveling component. Below 30mb the model amplitudes 
are In good agreement for wave 1 (compared to van Loon and Jenne, 1972), while 
waves 2 and 3 may be up to a factor of two too large In the mid-troposphere 



(van Loon and Jenne, 1972; Trenberth, 1980). The phase of wave 1 appears to be 
40® too far west, while the phase of wave 3 appears 40® too far east. In the 
Northern Hemisphere the amplitudes agree well with observations (van Loon et 
al., 1973; Defant et al., 1979; Reiter and Westhoff, 1981) with small 
amplitudes of about 20m for waves 2 and 3 in the mid-troposphere, and 
amplitudes for wave 1 of about 100m in the upper troposphere. The phases show 
strong variations with latitude and altitude, a commonplace occurrence for 
small amplitude waves, and evident in observations (Defant et al., 1979). In 
the stratosphere, the amplitudes remain quite small in the east wind regime. 
Note that the winter-summer difference in amplitude of the long waves at 500mb 
is substantially less in the Southern Hemisphere, in agreement with the 
observations of Trenberth (1980). 

t 

To summarize these comparisons: in both winter hemispheres wave 1 

appears to be well simulated in the troposphere and lower stratosphere, while 
waves 2 and 3 show some deviation from observations in the mid-troposphere. 
Given that the wave 1 forcing appears realistic, it should be possible to 
evaluate why both the model and observations show smaller wave 1 amplitudes in 
the Southern Hemisphere mid and upper stratosphere than is found in the 
Northern Hemisphere during winter. As shown in Fig. 20 and 21 there is little 
difference at 100 mb. Above lOOmb, however, wave 1 propagates better in the 
Northern Hemisphere. A complete discussion of wave propagation will be given 
in Section 3e. 

KELVIN WAVES (?): One other aspect of the wave spectrum relates to the 

discussion in Section 3b, subsection 2, and Fig. 10 which shows the semiannual 
zonal wind variation between 3.4 and 0. 17mb at the equator. It has been 
hypothesized (Hlrota, 1980) that the westerly phase of this variation is 
caused by the upward propagation of Kelvin waves, which provide westerly 
momentum when they are damped in the tropical stratosphere. Kelvin waves 
travel eastward, have long wavelengths, with zonal wind but no meridional wind 
oscillations. They should propagate best into the upper stratosphere when the 
zonal winds are easterly (during the solstices), and thus provide the momentum 
source for the change to west winds (during the equinoxes). 

Fig. 22 shows the time power spectrum of the zonal wind for eastward 
traveling components of wave number 1 near the equator averaged for the months 
around the solstices and equinoxes when the model tropical stratosphere winds 
were most easterly (July, August, January and February) and most westerly 
(April, May, September, October). The wave number one component was 
determined by Fourier analysis, and the time power spectrum was obtained using 
the maximum entropy method (e.g. , Hayashl, 1977) with order 10 regression 
(which was found to give the best results for artificially generated waves 
during a 30 day period). The results are shown for the upper troposphere and 
several levels in the stratosphere. As the winds did not obtain the proper 
west wind strength during the equinoxes it might be presumed that the model 
results will not be as informative as would otherwise have been possible, but 
certain similarities between the model simulations and theory are apparent. 

At 300mb and 50mb the energy at the solstices is only very slightly greater 
than that at the equinoxes; the east winds in the lower stratosphere in the 
tropics are also only slightly stronger during the solstices (Fig. 9b). By 
lOmb there is significantly more energy during the solstices, and this holds 
true at Imb. As altitude increases the power shifts to shorter periods, 
peaking in the 10-15 day period range by lOmb, with significant energy at 5-10 
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Fl^. 22. Time power spectrum (per unit mass) of the zonal wind for eastward 

traveling; components of wave number 1 near the equator. Results for 
the solstice (solid line) are for the months of January, February, 
July and August. Results for the equinox, dashed line, are for 
April, May, September and October. 


days by Imb. This is expected as the damping rate of Newtonian cooling 
Increases with decreasing Doppler-shifted frequency (Holton, 1975). 

Observations show Kelvin waves with periods of 10-20 days in the lower 
stratosphere (e.g. , Wallace and Kousky, 1968) and periods of 4-9 days in the 
stratopause region (Hlrota, 1980). It can also be seen that the longest 
periods occur at the solstices, so the Doppler-shifted frequency has less 
variation from season to season (as observed by Hlrota, 1980). 

Besides the above mentioned characteristics, these waves have the 
necessary attribute that the meridional velocity fluctuations are an order of 
magnitude less than the zonal. To truly Identify them as Kelvin waves 
requires analysis beyond the scope of this paper. The results shown provide 
"prima facie" evidence that Kelvin waves are propagating vertically and 
providing for angular momentum convergence (which is occurring due to vertical 
transports by eddies) In the model upper stratosphere. Transports will be 
discussed In Section 3e, as will the Influence on the equatorial mean flow of 
Rossby wave propagation from mid-latitudes, which Is thought to be responsible 
for generating the east wind phase of the semiannual soclllatlon. 

e. Wave propagation and transports 

1) HORIZONTAL EDDY TRANSPORTS 

SENSIBLE HEAT: Fig. 23 shows the northward transport of sensible heat 

by eddies during January and July. The tropospheric values In the Northern 
Hemisphere during January agree well with observations (Oort and Rasmusson, 

1971; Lau, 1979). (For further comparison with observations in the 
troposphere, see I). In the lower stratosphere, the values range from 
15K m s“l at lOOmb to 45K m s“^ at 30mb. This Is In agreement with the 
observations for January In this height range of Oort and Rasmusson (1971), 

Newell et al. (1969), Doppllck (1971), Lau (1979) and Crane et al. (1980); 
values shown by Geller et al. (1983) are almost a factor of two higher at 
30mb. This Is an Important diagnostic, as the cold lower stratospheric 
temperatures In the model might be caused by a deficiency In sensible heat 
transport, in addition to the radiative Impact of too much water vapor. 

Labltzke and Goretzki (1982) show that there Is great variability from year to 
year at 30mb. The underestimate of wave number 3 at lOOmb should not strongly 
Impact the sensible heat transport, if the analysis of Kao and Sagendorff (1980) 
for the winter of 1964 la representative. 

By lOmb the value has risen to 130K m s~^ in the model, in agreement with 
the observations of Geller et al. (1983) and slightly less than those reported by 
Doppllck (1971). The values then rise to a peak of almost 350K m s“^ at Imb, 
about 50% greater than that reported by Geller et al. (1983). This excess 
sensible heat transport Is consistent with the larger amplitude wave 1 produced 
by the model at that altitude. 

In the Southern Hemisphere troposphere during this month the model values 
are slightly too small to appear In a contour, but they are in agreement with the 
value of lOK m s“^ reported by van Loon (1980). Values remain small In the 
Southern Hemisphere summer stratosphere, as the eddy energy Is very small. 

In July, the small Northern Hemisphere sensible heat transports in the 
troposphere are in agreement with the observations of Oort and Rasmusson (1971) 
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of 5 k m s”^, while the small stratosphere values are again indicative of small 
summertime eddy energy. 

In the Southern Hemisphere winter troposphere the values are in agreement 
with results shown by van Loon (1980). The model values rise to a peak near Imb, 
with a value almost three times larger than that shown by Hartmann (1976a) for 
standing eddies. As noted earlier, the energy in the Southern Hemisphere winter 
stratosphere in the model is made up of 25% standing eddies and 75% transient 
eddies; the model transport associated with standing eddies is a little less than 
one-third the value shown in the figure, and thus in agreement with Hartmann. 

Thus the model's total poleward transport of sensible heat in the Southern 
Hemisphere winter stratopause region is dominated by the transient fluxes. 

WESTERLY hfiDMENTUM: The northward transport of westerly momentum by 

eddies for January and July is shown in Fig. 24. During January in the 
Northern Hemisphere troposphere the model simulation is in good agreement with 
the observations of Oort and Rasmusson (1971) and Lau (1979). The values in 
the low and mid stratosphere are in agreement with observations reported by 
Crane et al. (1980) although they may be slightly low compared with the 
observations of Geller et al. (1983); as shown by Labitzke and Goretzkl (1982) 
there is much year to year (and week to week) variation for January. It is 
not obvious from a comparison with observations that the model is 
underestimating the northward momentum transport in the region around lOOmb, 
which would affect the zonal wind distribution above the tropospheric jet. In 
the upper stratosphere the model is in fairly good agreement (within about 
20%) of the value reported by Geller et al. (1983) but it climbs to larger 
values, by almost a factor of two, in the lower mesosphere. This again is 
apparently due to the overestimate of the wave amplitudes at this height. 

In the Southern Hemisphere troposphere, there is only a small variation 
in the momentum transport from summer to winter, in both the model and 
observations (Oort, 1982). At higher levels during July, the Southern 
Hemisphere transports can be compared with observations of standing eddy 
transports shown by Hartmann (1976a). Once again the standing eddy transports 
are in agreement with the observations, while representing only one-third of 
the total eddy momentum transport. 

The July transports in the Northern Hemisphere are small, and may 
underest Imate^he transport in the upper troposphere by a factor of two 
(compared with Oort and Rasmusson, 1971; Oort, 1982). Note that the Northern 
Hemisphere tropospheric jet stream in July (Fig. 7) is also somewhat too 
weak. 


2) VERTICAL PROPAGATION AND TRANSPORTS 

GEOPOTENTIAL ENERGY: Shown in Fig. 25 is the vertical transport of 

geopotential energy by eddies for January and July. In January in both 
hemispheres the transport is downward from the middle troposphere to the 
Surface, and upward from the mid-troposphere into the stratosphere. It thus 
indicates that the energy generated in the middle troposphere due to 
barocllnlc conversions (see Fig. 17) is exported both upwards and downwards. 
The convergence of this flux is also shown in Fig. 17; this flux is often 
called the eddy pressure interaction term, or simply the boundary flux. The 
directions of transport are consistent with the results shown by Oort and 
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Rasmussen (1971) for standing eddies* with the implied vertical transport 
shown by Rung (1966), and in essentially complete agreement in both magnitude 
and direction with the transport shown by Manabe and Mahlman (1976) for their 
model January. 

The vertical transport through the lOOmb level in the Northern Hemisphere 
is similar to that shown by Manabe and Mahlman (1976) for the total eddy 
transport in their model; it Is 50% greater than the standing eddy flux 
computed by Newell and Richards (1969) and Oort and Rasmussen (1971), and 2 to 
3 times greater than the standing eddy flux obtained by Doppllck (1971), The 
standing wave energy at this altitude In the model is one-third the total eddy 
energy, so the total eddy vertical transport is not unreasonable. Between 
lOOmb and lOmb the model values are 30% excessive compared to the results 
shown by McNulty (1976) for standing wave fluxes. As shown in the figure 
there is only a gradual decrease of the vertical flux through the stratosphere 
indicating that waves are propagating relatively freely (especially wave 1) in 
the Northern Hemisphere. In the Southern Hemisphere there is a strong 
decrease in geopotential flux above 200mb, as will be discussed in greater 
detail below. 

During July, in the Northern Hemisphere the vertical transport above 
200mb decreases to very small values, in agreement with the relative change 
shown by Oort and Rasmussen (1971) for standing eddy transports. The vertical 
flux remains small throughout the summer hemisphere. In the Southern 
Hemisphere the vertical flux Is about a factor of two larger, from all eddies, 
than the values shown by Hartmann (1976a) for standing eddies; the Importance 
of transient eddies in the model's Southern Hemisphere stratosphere has 
already been mentioned* 

Comparing the two winter hemispheres, it can be seen that the eddy 
vertical flux is similar at 21mb, but decreases above in the Southern 
Hemisphere relative to the Northern Hemisphere. This is shown in another 
format below. 


ELIASSEN-PALM FLUX; Fig. 26 shows the Eliassen-Palm (EP) flux 
(Ellassen-Palm, 1961) for January and July, as well as the flux divergence. 

For presentation purposes constant length arrows are used in the troposphere; 
in addition, above lOOmb the vertical component is multiplied by a factor of 
250 relative to the horizontal, below lOOmb by a factor of 7. The dominant 
terms in this flux are the eddy transport of sensible heat, which provides for 
an upward EP flux when the transport is poleward, and the eddy transport of 
westerly momentum, which provides an equatorward EP flux when the transport is 
northward* (The computation which produced these figures includes the 
additional, non-quaslgeostrophlc terms as well, a heat transport term 
multiplying the shear in the northward component, and the vertical eddy 
momentum flux in the vertical component). For steady waves under adiabatic 
and frictionless conditions these fluxes are equivalent to the northward and 
vertical geopotential energy fluxes normalized by the zonal wind, and thus 
represent the direction of the group velocity for propagating waves in 
regions where the group velocity is well defined (Palmer, 1982). It is 
immediately obvious from this figure, compared with Fig. 25, that the vertical 
EP flux, which is dominated by the poleward sensible heat transport shown in 
Fig. 23, is opposite in direction to the actual geopotential flux in the low and 



Pressure (mb) Pressure (mb) 











55 


middle troposphere. A similar result can be seen in observations of stationary 
eddy transports (Oort and Sasmusson, 1971). In this region the assumptions of 
steady flow and adiabatic frictionless conditions are strongly violated; thus 
also Is the relationship between poleward heat transport and vertical 
geopotential transport discussed, for example, by Hartmann (1976b). (This is a 
similar result to that given In Section 3d, subsection 2, for the relationship 
between northward heat transport and baroclinlc eddy generation). 

From the upper troposphere through the stratosphere the flux Is aligned 
with the flux of geopotential energy (normalized by the zonal wind). It 
graphically indicates vertical propagation at mid and high latitudes, of greater 
magnitude In the winter hemispheres, consistent with the sensible heat transport 
shown In Fig. 23. In the winter hemispheres the flux Is equatorward from 
mid-latitudes, with a tendency for poleward flux at the highest latitudes - this 
is consistent «d.th the momentum transport relationship discussed above and shown 
In Fig. 24. In the low and mid-stratosphere the flux Is from the winter 
hemisphere Into the summer hemisphere, while the reverse appears true above 
about 2mb. 

The flux directions can be compared with those shown by Edmon et al. 

(1980) and Raroly (1982) for all eddies In the troposphere, with observations 
of Sato (1980) and Geller et al. (1983) for troposphere and stratosphere, and 
with those of Smith (1983) In the stratosphere for wave 1. In the troposphere 
the directions are In complete agreement except that the summer subtropical 
fluxes in the Northern Hemisphere are not equatorward In the model (remember 
that the westerly momentum transport In this region during this month, as shown 
In Fig. 24, is underestimated, as Is the strength of the summertime subtropical 
jet). The magnitudes are presumably realistic, (although it Is hard to tell 
from the diagrams In the literature), because as discussed above, the sensible 
heat and momentum transports are generally In agreement with observations from 
the troposphere through the middle stratosphere. 

REFRACTIVE INDEX: The characteristics of wave propagation are understood 

In the quaslgeostrophlc system In terms of the refractive Index squared (as 
defined by Matsuno, 1970, and modified by Palmer, 1982). Waves will not 
propagate through (infinitely thick) regions of negative refractive Index, 
which Is caused either by a negative latitudinal gradient of quasl-geostrophlc 

potential vortlcity ( qy), east winds or very strong west winds; the largest EP 
fluxes occur In regions of moderate positive values of the refractive Index, 
curving toward the direction of increasing Index In regions of large gradients 
(Smith, 1983). 

Fig. 27 shows the refractive Index-squared for wave number 1 for January 
and July, as modified by Palmer (values divided by sln^ 6 where 0 Is the 
latitude). The value for wave 1 Is shown since it provides the greatest 
portion of energy In the stratosphere; thus while an exact comparison with the 
total wave propagation shown by the EP flux diagrams Is not possible, a 
qualitative relationship can be determined. The refractive index-squared Is 
negative In the summer hemisphere stratospheres above about 30mb due to the 
east winds; the vertical EP flux (and the vertical transport of geopotential 
energy) decreases to very small values In the summer hemisphere above that 
altitude. At low latitudes the refractive Index-squared becomes very large and 
negative (due to east winds and the latitudinal weighting), or occasionally 
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Fig. 27. Latitude-height profile of the modified refractive Index-squared of 
wave number 1 for (a) January and (b) July. Regions with values 
greater than 100 are striped, with values <-100 are shaded. 





very large and positive so the EP flux is basically horizontal. The EP flux 
is large in regions of positive refractive index-squared, and is guided into 
the core of the west winds (as shown in Fig. 7) due to the relatively large 

value of qy there (which in turn is due to the curvature of the latitudinal 
wind profile). This is all in accordance with expectations and observations, 
including the negative value of the refractive index-squared found at the 
highest latitudes in winter above lOmb (Smith, 1983). 

However, there are some discrepancies: the diagnostic shows a region of 

negative refractive index-squared, due to a negative latitudinal gradient of 
quasi-geostrophlc potential vortlclty, in the mid to upper troposphere at 
upper mid-latitudes of the winter hemisphere. Some observations show negative 

values of qy in the mid-troposphere at these latitudes (Fullmer, 1982), 
although others give only positive values of the refraction index-squared 
(Lin, 1982; Karoly, 1982). As shown by the vertical EP flux, the waves do not 
appear to be making any special effort to avoid the region. Furthermore, 
observations (Lin, 1982) imply a minimum in the refractive index-squared 
directly above the tropospheric jet, which does not occur in the model, 
although one does occur around 10 mb. These negative values of qy are caused 
by the failure of the model to decrease the latitudinal temperature gradient 
with increasing altitude in the lower polar stratosphere. This deficiency is 
equivalent to the inability of the model to provide a decrease in the zonal 
wind at mid-latitudes in the lower stratosphere, as discussed in Section 3b. 

As shown by Fullmer (1982), qy is very sensitive to changes in the zonal wind. 

It is curious, however, that the EP flux does not seem to be affected strongly 
by these deviations, %diich raises the question of the applicability of this 
quasi-geostrophlc diagnostic to the general circulation model, at least 
Insofar as details are concerned. The index is slightly larger in winter in 
the Northern Hemisphere above 20mb than in the Southern Hemisphere. As noted 
above, there is considerably less wave propagation in the Southern Hemisphere 
at these heights. The zonal winds are some 20 ms~^ stronger at these 
altitudes during winter in the Southern Hemisphere and this would act to limit 
vertical propagation all other things being equal. 

DIVERGENCE OF EP FLUX: Also shown in Fig. 26 is the divergence of the 

EP flux; this is a measure of the eddy effect on the mean flow, being zero for 

steady waves Ih adiabatic frictionless flow with no critical level (Andrews 
and McIntyre, 1976). In the upper troposphere, the zero wind line, a critical 
line for stationary waves, occurs at low altitudes of the winter hemisphere; 
the EP convergence there acts to decelerate the zonal flow. The same effect 
is noted in the stratopause region, and this, at least in the model, verifies 
the contention of Dickinson (1969) and Hopkins (1975) that the east wind phase 
of the semiannual oscillation appears to be caused 1^ absorption of planetary 
waves propagating from the winter hemisphere. Note that in Fig. 24 there is 
strong westerly momentum divergence in the subtropics of the winter hemisphere 
due to eddy transports. The large values of the EP flux divergence near the 
stratopause in the region of relatively large refractive index-squared may be 
because the group velocity is decreased in this region and the vertical scale 
of 'the waves decreases (relative to regions of moderate refractive 
index-squared), thus dissipation is larger (Smith, 1983). The convergences in 
the upper levels of the winter hemispheres are due to the effect of the 



parameterized drag on the waves. In regions of EP flux divergence the waves 
are accelerating the mean flow - note the strong eddy Induced westerly momentum 
convergences shown In these regions In Fig. 24. The regions of divergence In 
the west wind maximum, the region of convergence in the subtropical winter 
stratopause region, and also the region of convergence at the same altitude In 
the summer hemisphere subtropics have all been Identified (Geller et al., 

1983). The tropospheric convergence region In the winter subtropics Is in 
agreement with observations (Karoly, 1982; Edmon et al., 1980; Geller et al. , 
1983) and the entire tropospheric pattern at mid-latitudes Is qualitatively 
similar to the picture given by Edmon et al. (1980) for the time-averaged life 
cycle of baroclinlc waves. 

In the summer hemisphere stratospheres there Is little In the way of EP 
flux divergence. This represents the fact that the eddies are not affecting 
the mean flow In these east wind regimes. 

Experiments calculating the transport of passive tracers in the 
stratosphere (e.g.. Hunt and Manabe, 1968) have Indicated that the net 
transport Is a small difference between transports In one direction due to 
eddies, and transports In the opposite direction by the mean meridional 
circulations. This Is now recognized as an example of a type of 
non-interaction theorem (Mahlman et al. , 1980), as the mean cells are also set 
up by the eddies. The cancellation will be Incomplete to the extent that 
conditions similar to those for the noninteraction theorem are violated, l.e.. 
In regions where the EP flux divergence Is nonzero. While this holds true for 
passive tracers. It also applies to atmospheric parameters. In the regions of 
the free atmosphere In the model (those not strongly Influenced by friction 
like that occurring In the lower troposphere or the upper three levels), 
winter hemispheric transports occur essentially In regions of nonzero EP flux 
divergence. In the region of EP flux divergence the total atmospheric 
horizontal transport (by eddies plus the mean circulation) of dry static 
energy (sensible heat plus geopotential energy), and of angular momentum Is 
equatorward; in regions of EP flux convergence, the fluxes are poleward. In 
the summer hemispheres there Is little divergence and little transport. 

POTENTIAL VORTICITY: One additional vertical transport which Is of 

interest Is that of potential vorticlty. Mixing between the troposphere and 
stratosphere demands that the higher values of potential vorticlty In the 
stratosphere be brought down Into the troposphere, and be dissipated by 
friction or nonadlabatlc effects. The scale at which this transport occurs is 
not well known, and there Is some question whether It can be adequately 
resolved in a coarse grid model. With potential vorticlty written in pressure 
coordinates (Gldel and Shapiro, 1980), the vertical transport through the 
lOOmb layer for January and July Is shown in Fig. 28, for both eddies and the 
total (eddies plus the mean circulation). The locations of downward 
transport, strongest In mid-latitudes of the winter hemisphere, but also 
occurring In upper mid-latitudes of the summer hemisphere agree with the NCAR 
GCM results reported by Gldel and Shapiro (1980), and, for the Southern 
Hemisphere winter, the analysis of Hartmann (1977). The eddy transports are 
generally the same sign as the total transports, but In winter are not the 
dominant component. The validity of using a coarse grid model to Investigate 
tropospheric-stratospheric transports can be tested by using the model to 
advect passive tracers and comparing with appropriate observations. 
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Fig. 28. Vertical transport of potential vortlclty through the lOOmb level by 
eddies and total for January (top) and July (bottom). 




f . Time sequences 


As noted previously, the Northern Hemisphere stratosphere underwent 
pulsing events during the winter, with warming amplitudes most noticeable 
during January. A sequence of events Is shown in Fig. 29 at 0.75mb (45-50km 
in the Northern Hemisphere) starting on January 1, the 14th month of the 
simulation. The left hand side of each row shows the geopotential height map 
averaged over the day, while the right hand side shows the temperature 
pattern. Every fourth day Is shown through day 13. Starting on January 1, 
the first maps Indicate that the polar cyclone was situated over Greenland, 
with a ridge centered over the Mediterranean Sea. The warmest temperatures, 
+5®C, were located in a warm region just to the northwest of the high, while 
the coldest temperatures were found In the Aleutian area. Four days later the 
ridge had strengthened by 1km, and was now centered somewhat further east, and 
slightly north, at SO^N, 90“E. The polar cyclone had moved slightly westward. 
The warmest region had strengthened and moved far to the north, with 10®C 
located north of the ridge near the pole. By day 9 the ridge had weakened and 
moved eastward, to a position over the mld-Paclflc, while the polar low had 
also moved eastward, to the side of the pole near the Greenwich meridian. 

Both the warm and cold regions showed little movement, although the warm 
region was weaker, now at 0°C. On day 13 the high Is now again further west, 
near 120 ®E, and somewhat stronger, while the low Is still east of Greenland. 
The warm region has moved south, and Intensified strongly, now at 20®C, while 
the cold region has moved eastward to a position over the Canadian Marltimes, 
and is 5®C colder than it had been. 

The movements described here In the upper stratosphere are typical of 
those seen at various levels In the stratosphere associated with major or 
minor warmings (e.g., Williams, 1968). There are several points to be made. 
There was no major warming, since the polar cyclone never retreated from Its 
near pole position, and the latitudinal temperature gradient did not reverse 
(and further the definition of a major warming requires that this happen at 
lOmb). While wave 1 dominates the maps, the event Is not restricted to 
standing or slowly moving wave 1 processes - note the warm and cold regions on 
day 13, at their greatest amplitudes, were separated by just 120® longitude at 
60®N. The whole pattern of shifting locations for the ridge and trough could 
be analyzed in terms of travelling and standing waves. The time power 
spectrum for geopotential at 50®N and Imb shows significant energy at 15-30 
day periods for westward traveling waves. Finally, if one observes the 
temperatures near the equator at the approximate longitude of the warm region. 
It can be seen that when the region at high latitudes is warmest, the 
equatorial region Is coldest. When the warm region Is at 5®C, the equatorial 
temperature to the south Is at -10®C (day 1); when It is at 10®C, the 
equatorial temperature Is at -15®C (day 5); and when it Is at 20®C, the 
equatorial temperature is at -20®C (day 13). This simultaneous negative 
correlation has been observed during warmings (e.g., Labltzke and Barnett, 
1973), and has been associated with wave transports set up as planetary waves 
are absorbed (Matsuno, 1971); the mechanical forcing sets up a residual mean 
circulation (Plumb, 1982). 

The appearance of these warming pulses as a function of time and altitude 
can be seen by looking at the zonal average temperature at 67 *N from late 
December through January. This is shown In Fig. 30. The amplitude of the 
perturbation Is reduced somewhat by using the zonal average, as can be seen In 
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Fig. 29a 


Northern Hemisphere geopotential heights (left) and temperatures 
(right) for the 0. 75mb level for (a) January 1; (b) January 5. 
Height contours are In meters above or below 45 km. 









OR/GIWAL PAGE fS 
or POOR QUALITY 



Days 


Fig. 30. 


Temperature deviation from the 35 day average at each pressure level 
as a function of day from December 27 through January 31 for the 
zonal mean temperature at 67 *N. 



Fig. 29. From the figure it appears as if the warming pulses originate at 
about 65km and propagate down to 45km in 4-5 days. Propagation to 25km occurs 
more quickly, in 1-3 days. In the stratosphere the perturbation reaches its 
largest amplitude around 42km and is hard to detect at levels below 20km. 

This type of downward propagation of warming events is a well-known 
phenomenon. It is possibly associated with vertically propagating transient 
waves initiating their effect at high altitudes, where the density is low and 
wave amplitude large, and altering the atmosphere so the effect is felt 
subsequently at successively lower levels (with or without the formation of 
critical levels -e.g., Matsuno, 1971), Note also that in association with the 
phase lag, the warming in the upper stratosphere occurs simultaneously with 
cooling in the mesosphere, also an observed feature (e.g., Scott, 1972) 
in agreement with Matsuno 's model and the properties of the forced Hough 
functions (Plumb, 1982). The pulses in the model occur with a ~ 

quasi-periodicity of 9-13 days, in agreement with observations of 1-3 week 
pulsations (Madden, 1975) (and the time power spectrum at 50“N), and may be 
thought of as indicating the presence of traveling wave modes which 
alternately constructively and destructively interfere with the 
quasl-stationary waves (Madden, 1975). 

It is beyond the scope of this paper to analyze these events in great 
detail. The time sequence is presented here to indicate the model's fidelity 
to observations. It is important for modeling changes in the stratosphere 
associated with climate perturbations that the model be able to reproduce 
transient features as well as climatological averages, for they may lead to 
major or minor warmings which impact the transport of constituents such as 
ozone. 

4. Discussion 

The comparisons with observations shown in the previous section indicate 
that the model is capable of generating a realistic troposphere and 
stratosphere. The model produces generally realistic temperatures, heights, 
winds, long waves, energy fluxes and transports, as well as realistic 
latitudinal, longitudinal, seasonal and hemispheric variations of the above 
quantities. The comparison also shows the need for improvements in the 
following areas: 

(1) The temperature in the Northern Hemisphere winter polar lower 
stratosphere is too cold by some 5-10*C. This appears to be mostly the result 
of the excessive water vapor which is allowed to accumulate in the region and 
increase long wave radiation losses. The obvious solution is to allow water 
vapor to condense above the lOOmb level, a process which was not included in 
the current version of the model. 

(2) As a result of these cold polar temperatures, or at least consistent 
with them, the winds in the lower stratosphere are too strong. The resulting 
latitudinal wind profile helps to produce an apparently unrealistic refraction 
index profile, although the effect of this on the long wave propagation is not 
overly apparent. If the problem results from the excessive water vapor 
cooling it should be alleviated by the correction suggested above. 

(3) The apparent deficiency in wave 3 energy in the troposphere may be 
related to the unrealistic refractive index profile, as suggested by Tenenbaum 



(1982). The magnitude of the underestimation is uncertain, however, since the 
alternate January has significantly more wave 3 energy. Wave 4 energy also 
appears to be deficient in the troposphere. 

(4) Temperatures at 50km and above in equatorial regions are too cold. 
This is probably due to the Inability of the current radiation scheme to 
properly simulate CO 2 long wave cooling above about 45km. The radiation 
scheme needs to be altered to remove this possible deficiency. 

(5) The excessively cold tropical stratopause temperatures result in a 
low latitude temperature gradient which is too large, and, from the thermal 
wind relationship, this results in the model's simulation of east winds being 
about two times too strong in the tropical lower mesosphere. 


(6) There is some evidence that the drag incorporated in the top three 
model levels is too strong in winter and could be reduced, so that the winter 
west winds would not decrease as rapidly with altitude above the stratospheric 
jet. The summer mesosphere winds, however, do not appear excessive. 

(7) Planetary wave amplitudes in the upper stratosphere and lower 
mesosphere are too large. The model does not Include photochemical damping 
associated with the temperature dependence of ozone reactions, which would 
likely reduce the amplitudes* A parameterized form of photochemical damping, 
(e.g.. Pels et al., 1980) should be included. 


(8) The possibility of Including explicit diffusion should be reexamined. 
The model currently uses none; obviously the choice of appropriate diffusion 
coefficients could by Itself alleviate some of the deficiencies noted above 
(for example the large planetary wave amplitudes in the upper stratosphere). 
The problem with using diffusion is that the uncertainty in the true magnitude 
of the diffusion coefficients provides leeway for using values which "tune” 
the model to the current climate, and thus limit the nature of the 
stratospheric reaction to climate changes. 

(9) The model produces apparent inertial oscillations in the tropical 
upper stratosphere. While these may exist, the model simulation may be 
resolution dependent, with a coarse resolution magnifying the effects of the 
instability. 

(10) This last point leads to the more general conclusion that model 
results may be affected by Increasing the vertical and horizontal resolutions. 
The modeled processes in the lower stratosphere, and 

tropospheric-stratospheric exchange may well differ with higher resolution, 
and it is important to quantify these effects. The coarse resolution model 
employed here is useful for testing and model development, and to the extent 
that it simulates climatology it is an adequate model for understanding the 
current stratosphere, but its use in climate change experiments must 
ultimately be validated against higher resolution versions. 

(11) Finally, the model needs to be run for at least several years to 
determine its true climatology and Interannual variability. 
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